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INTRODUCTION 


The  adverse  effects  of  changing  temperatures  on  electrical 
network  performance  is  of  increasing  importance.   Circuits  of 
high  linearity  and  stability  are  desired,  yet  they  must  often 
operate  under  extreme  temperature  variations.   Designers  may 
attempt  to  limit  adverse  temperature  effects  in  one.  of  two  ways. 
First,  the  circuit  may  be  designed  in  such  a  way  that  variations 
in  temperature  have  as  little  effect  as  possible.   This  is  usually 
only  a  partial  solution  and  is  costly  in  terms  of  desired  circuit 
performance.   Second,  additional  circuitry  may  be  designed  to  com- 
pensate for  the  undesired  temperature  effects.   For  example  the 
gain  of  a  feedback  amplifier  may  vary  with  temperature.   This 
variation  may  be  partially  eliminated  by  increasing  the  feedback, 
but  the  gain  is  decreased  in  the  process.   The  temperature  in- 
duced gain  variation  may  be  eliminated  by  making  the  feedback 
voltage  vary  with  temperature  in  the  proper  manner.  The  frequency- 
temperature  characteristic  of  a  crystal  resonator  may  be  linear- 
ized by  varying  the  crystal  load  capacitance.   This  capacitance 
may  be  provided  by  a  varactor  biased  by  the  proper  d.c.  control 
voltage.   In  both  of  these  examples  as  well  as  in  other  compensa- 
tion problems,  generation  of  a  temperature  variable  control  vol- 
tage is  required. 

The  subject  of  this  study  is  the  design  of  a  circuit  to 
provide  a  specified  temperature  compensated  control  voltage. 
The  possibility  of  approximating  a  thermistor  by  the  first  few 


terms  of  a  power  series  was  studied.   This  did  not  yield  useful 
results.   A  linear  approximation  was  made  for  a  thermistor  and 
resistor  in  parallel.   This  linear  approximation  was  used  in  the 
study  of  a  voltage  divider  and  a  bridge  network.   A  method  of 
generating  temperature  compensated  voltages  using  the  bridge  net- 
work was  devised.   This  generation  method  was  applied  to  the  pro- 
blem of  compensating  a  crystal  resonator. 


DESIGNING  A  TEMPERATURE  COMPENSATING  CIRCUIT 

The  design  of  a  temperature  compensation  circuit  requires 
a  device  whose  physical  properties  vary  with  temperature.   De- 
vices with  temperature  variable  resistance,  capacitance,  or  in- 
ductance are  available  or  could  be  constructed.   Certain  proper- 
ties are  necessary  or  desirable  in  such  a  device.   The  thermistor 
satisfies  most  of  them.   The  range  of  resistance  variation  with 
temperature  is  large.   Thermistor  resistance  can  vary  by  a  factor 
of  10   over  a  200°  C  temperature  range.   The  thermistor  is  also 
inexpensive.   The  effect  of  the  thermistor  on  the  compensated 
circuit  or  device  can  be  easily  modified  by  placing  it  in  various 
resistance  networks.   If  a  device  were  compensated  by  means  of  a 
temperature  variable  capacitance  or  inductance,  the  effect  of  the 
compensating  element  would  have  to  be  modified  by  capacitive  or 
inductive  networks.   The  resistive  network  necessary  for  a  ther- 
mistor would  be  less  expensive  than  an  inductive  or  capacitive 
network. 

Thermistor  Characteristics 

Thermistors  are  semiconductor  devices.   Their  resistivity 
lies  between  that  of  conductors  and  insulators.   Thermistors  are 
made  from  both  N-type  and  P-type  semiconductor  materials, 
generally  mixtures  of  manganese,  nickel,  and  cobalt  oxides.   Con- 
duction in  either  type  of  material  is  the  result  of  the  generation 
of  intrinsic  electron-hole  pairs.   These  are  formed  by  thermally 


excited  electrons  being  elevated  to  the  conduction  band.   At 
higher  temperatures  more  electron-hole  pairs  are  formed  and  hence 
there  is  higher  conduction  and  lower  resistance.   The  resistance 
of  a  thermistor  is  given  by 


R   -  R   e 

t     o 


BC1/,,,  -  X/T  ) 


(1) 


where  R   is  the  resistance  at  some  nominal  temperature  I  ,  T  is 
o  C 

the  temperature  of  interest,  and  B  is  a  constant  for  a  particular 

thermistor.   Depending  on  thermistor  materials  and  construction 

R   may  vary  from  a  low  of  100  ohms  or  less  to  a  high  of  several 
o 

megohms.   Not  strictly  a  constant,  B  may  vary  slightly  with 
temperature  but  may  be  considered  constant  over  any  150°  C  temp- 
erature range.   Thermistor  materials  display  a  wide  range  of  B 
values,  varying  from  a  low  of  almost  zero  to  a  high  of  about 
7,000.   The  units  of  B  are  degrees  centigrade  or  Kelvin.   This 
wide  range  of  thermistor  constants  makes  possible  the  selection 
of  a  thermistor  to  suit  the  needs  of  a  particular  compensation 

problem.   The  thermistor  temperature  coefficient  of  resistance, 

dR  /dT 
which  may  be  derived  from  Eq.  (1)  by  evaluating  — g ,  is 


V     • 


(2) 


The  temperature  coefficient  of  resistance  is  negative  as  would  be 
expected  of  a  semiconductor  device  and  is  highly  nonlinear  with 
respect  to  temperature. 
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When  placed  in  a  practical  circuit,  a  thermistor  will  have 
some  voltage  across  it.   A  current  will  flow  causing  power  to  be 
dissipated  in  the  thermistor.   This  power  dissipation  will  result 
in  an  increase  in  the  thermistor  temperature  and  consequently  a 
decrease  in  resistance.   A  logarithmic  plot  of  a  static  voltage- 
current  curve  is  given  in  Fig.  1.   The  numbers  along  the  curve 
give  the  rise  in  temperature  above  ambient.   A  current  of  less 
than  0.2  milliamperes ,  resulting  in  less  than  a  2  degree  tempera- 
ture rise  has  little  effect  on  the  resistance  of  this  thermistor. 
The  dissipation  constant  of  a  thermistor,  C,  is 


/ 


AT 


(3) 


where  P  is  the  number  of  watts  of  power  dissipated  by  the  ther- 
mistor and  AT  is  the  temperature  rise  above  ambient.   A  quantita- 
tive measure  of  the  effect  of  heating  on  thermistor  resistance  is 
the  power  sensitivity.   The  power  sensitivity  is  the  number  of 
watts  of  power  which  may  be  dissipated  before  the  thermistor 
resistance  decreases  by  one  percent.   The  power  sensitivity,  S, 
is  expressed  in  terms  of  the  power  dissipation  constant  and  the 
temperature  coefficient  of  resistance  by 


a  x  100 


(4) 


The  dissipation  constant  is  affected  by  the  surface  area  and  the 
volume  of  the  thermistor,  the  nature  of  the  surface,  and  the 
thermal  conductivity  of  the  surroundings  and  supports.   The 


dissipation  constant  may  be  altered  by  varying  the  physical  size 
and  shape  of  the  thermistor,  however  this  may  also  affect  the  value 

of  R  .   The  time  necessary  for  the  thermistor  to  respond  to  a 

o  '  * 

change  in  ambient  temperature  is  inversely  proportional  to  the 
dissipation  constant.   An  average  value  for  C  is  5  milliwatts  per 
degree  centigrade,  however  this  is  variable  by  a  factor  of  10  or 
greater . 

When  used  for  a  temperature  compensation  element,  the  tempera- 
ture and  hence  the  resistance  of  the  thermistor  must  be  deter- 
mined by  ambient  conditions  and  not  internal  heating.   A  necessary 
restriction  for  such  compensation  circuits  is  that  thermistor 
currents  be  limited  to  low  levels.   Values  of  currents  well  within 
those  dictated  by  the  rated  power  sensitivity  should  be  used. 
Since  ambient  conditions  may  change  rapidly,  the  dissipation 
constant  should  be  chosen  so  that  the  compensating  thermistors 
and  the  compensated  elements  both  respond  to  a  change  in  ambient 
conditions  at  approximately  the  same  rate.   If  this  is  not  done, 
sudden  changes  in  ambient  conditions  will  cause  transient  periods 

of  under  or  over  compensation.   Since  both  C  and  R   are  a  function 

o 

of  thermistor  size,  it  might  be  impossible  to  satisfy  specified 
restrictions  on  both  parameters  in  a  single  thermistor.   This 
problem  might  be  eliminated  by  using  the  desired  value  of  C  to 
determine  the  thermistor  type  and  by  placing  two  or  more  units 
in  series  or  parallel  to  obtain  the  desired  equivalent  value  of 
R  .   Finally,  it  is  important  that  the  compensated  elements  and 
the  compensating  thermistors  be  placed  in  close  physical  proximity, 
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since  the  temperature  difference  possible  between  two  network 
points  is  not  necessarily  constant  in  time.   For  example,  a 
sudden  temperature  increase  in  the  vicinity  of  a  certain  resistor 
caused  by  an  increase  in  load  current  would  not  immediately  re- 
sult in  a  temperature  increase  at  a  more  distant  point  in  the 
network . 

Several  temperature  compensation  methods  require  the  use  of 
a  d.c.  control  voltage  whose  amplitude  varies  with  temperature. 
Such  a  voltage  is  illustrated  in  Fig.  2,  where  temperatures  T   and 

T.  are  the  end  points  of  the  range  of  desired  compensation.   Often 
b 

this  voltage  is  given  or  may  be  accurately  approximated  over  a 
desired  temperature  range  by  a  polynomial  f(T).   Such  a  control 
voltage  may  be  realized  by  a  resistive  network  having  a  tempera- 
ture variable  transfer  function,  K(T),  driven  by  a  constant  vol- 
tage source.   The  transfer  function  may  in  general  be  approximated 
by  a  ratio  of  polynomials  such  that, 


K(T) 


N(T) 
D(T) 


(5) 


This  network  would  include  thermistors  as  the  temperature  varia- 
ble elements  and  would  hopefully  be  such  that  K(T)xV,  V  being  the 
source  voltage,  is  a  close  approximation  of  f(T)  over  the  limited 
range  of  interest,  T   to  T,  . 
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Power  Series  Representation  of  a  Thermistor 


Equation  (1)  may  be  written  as  a  power  series  in  T  expanded 

about  some  temperature  T  »  q.   It  is  quite  possible  that  the  first 

few  terms  of  the  power  series  might  give  an  accurate  approximation 

of  R   over  the  temperature  range  T   to  T,  .   The  power  series 
t  a      b 

representation  of  the  thermistor  might  then  be  used  in  the  trans- 
fer function  equations  of  a  network.   This  transfer  function  is 
then  of  the  form  K(T) .   The  known  coefficients  of  f(T)  and  the 
coefficients  of  the  power  series  provide  sufficient  information 
to  determine  the  unknown  coefficients  of  K(T).   The  coefficients 
of  K(T)  would  then  be  used  to  find  the  values  of  the  network 
resistors  which  are  not  temperature  variant. 

To  determine  the  accuracy  of  the  power  series  representation 
of  R  ,  the  power  series  was  computed  using  the  first  seven  terms. 
Restating  Eq.  (1), 


-B.     B/        B.j 

R   =  R   e     oe     =ke 
t     o 


(6) 


The  power  series  expansion  at  T 


q  is 
,2 


d,  x  (T-q)    d   x  (T-q) 

Rt  -  do +  — n +  — 2i + 


d,  x  (T-q)" 
+  -1— 61 + 


(7) 


Each  d.  of  Eq.  (7)  denotes  the  i    derivative  of  Eq.  (6)  evaluated 
at  T  =  q..   As  listed  below  the  first  seven  are 
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The    power    series    for    R      may    therefore    be    written    as 

2         j3         j4         j5         ^o5 
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An  iterative  computer  program  was  written  to  calculate  the 
thermistor  resistance  for  a  given  temperature  range  by  means  of 
the  power  series,  Eq.  (9),  and  the  defining  equation,  Eq.  (1). 
This  program  and  the  results  of  these  calculations  are  compiled 
in  Appendix  A.   The  data  was  computed  for  a  100°  K  temperature  range 
centered  at  q  "  283°  K.   A  comparison  of  the  two  sets  of  data 
showed  that  the  resistance  as  computed  by  the  power  series  was 
about  5  percent  low  at  233°  K  and  about  100  percent  high  at  333° 
K.   An  error  of  this  magnitude  made  the  power  series  approach 
unusable.   The  accuracy  of  the  power  series  could  be  increased 
by  adding  more  terms  to  the  series,  however,  this  would  have  made 
its  inclusion  in  the  network  equations  impractical. 

Linearized  Thermistor  Network 


A  linearized  temperature  variable  element  having  the  form 
R   =  m  -  nT,  where  m  and  n  are  positive  real  constants,  would  be 
desirable.   A  linearized  thermistor  element  would  lend  itself  to 


13 


intuitive  reasoning  in  the  design  of  a  practical  network  and  would 
greatly  simplify  the  network  equations.   A  clue  to  a  linearization 
method  for  thermistors  was  obtained  from  Becker,  et.  al.,  (1) 
in  their  discussion  of  temperature  compensation  of  copper  devices. 
If  a  thermistor  and  a  fixed  resistor  are  paralleled,  the  resistance 
of  the  combination  at  low  temperatures  is  limited  by  the  fixed 
resistor  and  at  high  temperatures  by  the  thermistor  resistance, 
which  is  close  to  zero.   This  is  shown  by  Fig.  3.   The  curve  of 
Fig.  3  is  essentially  linear  over  some  temperature  range  A.   By 
a  proper  choice  of  the  thermistor  and  resistor  combination  the 
linear  region  might  be  made  large  enough  to  span  the  desired  com- 
pensation range.   The  linearized  combination  could  then  be  used 
in  the  design  of  a  compensating  circuit  so  long  as  the  designer 
remembered  that  the  linear  approximation  is  valid  only  over  a 
limited  temperature  range.   A  trial  and  error  search  was  insti- 
tuted using  the  computer  program  of  Appendix  B  and  the  plotting 
routine  of  Appendix  C.   A  parallel  resistor  of  value  equal  to  the 
R   value  of  the  thermistor  gives  the  largest  linear  region.   As 
illustrated  by  Fig.  4,  the  slope  of  the  linear  region  of  the  curve 
can  be  varied  by  using  different  values  of  the  constant  B.   The 
linear  temperature  range  is  smaller  for  larger  values  of  B.   A 
piecewise  linear  approximation  can  be  made  for  the  whole  tempera- 
ture range.   Each  of  the  three  major  regions  can  be  approximated 
by  a  linear  section  as  shown  in  Fig.  5.   Regions  1  and  3  are  use- 
ful for  determining  marginal  behavior  of  a  network  and  for  deter- 
mining trial  values  after  a  network  has  been  selected.   Region  2, 
approximated  by  R   =  m  -  nT,  is  useful  in  the  design  of  a 
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compensation  network.   The  approximation  of  region  2  is  accurate 
within  about  +  5  percent,  while  the  high  and  low  temperature 
regions  are  accurate  within  only  about  +  20  percent. 

As  stated  previously,  a  voltage  of  the  form  of  Fig.  2  may 
be  obtained  by  operating  on  a  fixed  source  voltage  with  a  2-port 
network  whose  transfer  function  is  temperature  dependent.   In 
order  to  cope  with  many  practical  compensation  problems  the  out- 
put voltage  should  approximate  a  polynomial  of  at  least  third 
order  in  T  with  two  relative  extrema.   Two  operations,  amplifi- 
cation and  shifting  can  be  performed  on  the  output  voltage  to 
bring  it  into  line  with  the  desired  control  voltage.   Amplifica- 
tion is  in  the  normal  sense  the  multiplication  of  the  output 
voltage  by  a  constant  value.   Shifting  is  the  addition  of  a 
positive  or  negative  voltage  in  series  with  the  circuit  output 
or  amplified  output.   This  is  illustrated  by  the  diagram  and 
curves  of  Fig.  6.   In  Fig.  6b  the  curve  intersects  the  zero 
axis  at  three  locations.   These  points  are  known  as  zeros  of  the 
function.   The  number  of  zeros  and  their  location  determine  and 
help  to  describe  the  characteristics  of  the  curve.   As  shown  by 
Fig.  6c  shifting  is  a  versatile  tool  for  altering  a  voltage- 
temperature  curve.   Shifting  can  change  both  the  number  of  zeros 
and  their  locations.  Designing  a  network  which  will  have  one  zero 
and  a  relative  minimum  in  the  positive  region  as  in  Fig.  6c  is 
rather  difficult,  while  designing  a  region  as  in  Fig.  6b  is  much 
easier.   The  idea  of  shifting  makes  it  possible  to  design  a 
circuit  with  an  output  curve  like  Fig.  6b  and  easily  convert  it  to 
one  with  an  output  curve  like  Fig.  6c. 
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Fig.    6a.    Aroplifing    and    shifting    circuit    for    V 
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Fig.  6b.  Effects  of  amplification  on  the  voltage- 
temperature  plot. 
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Fig.  6c.  Effects  of  amplification  and  shifting  on  the 
voltage-temperature  plot. 
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Determining  a  Compensation  Network 

In  an  effort  to  find  a  network  with  sufficient  generality 
to  compensate  various  problems  with  only  a  change  of  circuit 
constants,  two  network  types  were  examined.   The  first  type 
relied  on  the  basic  voltage  divider  of  Fig.  7,  where  the  cir- 
cled resistors  with  enclosed  T  represent  linearized  thermistor 
elements.   The  linear  approximation  of  the  thermistor  was  used  in 
the  analysis  of  these  networks.   The  transfer  functionsof  these 
networks  are 


R+m-nT 


and 


m-nT 
R+m-nT 


(10) 


An  intuitive  view  of  these  transfer  functions  is  given  by  the 
voltage  transfer  ratio  vs.  temperature  plots  of  Fig.  8.   These 
results  seemed  to  lend  promise  of  obtaining  the  desired  transfer 
function  from  a  network  of  this  form.   The  next  step  was  the 
consideration  of  a  ladder  network  formed  by  cascading  these  vol- 
tage divider  sections.   This  network  is  shown  in  Fig.  9.   The 
transfer  function  of  the  network  of  Fig.  9  is  obtained  by  solving 
the  matrix  equations, 


(q+m-nT)     -(m-nT)  0 

-(m-nT)    (r+2m-2nT)     -(m-nT) 
0        -(m-nT)      (s+2m-2nT) 


(11) 
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Two  simple  voltage  divider  networks. 
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Voltage  transfer  functions  for  Fig.  7a  and  7b, 
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.  (13) 
Equation  (13)  is  of  the  general  form, 


e2T  +  ex 


3       2 
-c.T   +  e,T 
4       3 


(14) 


c2T  +  cx 


where  e   and  c   are  constants,  depending  on  the  values  of  m,  n, 
r,  q,  and  s.   Since  the  denominator  of  Eq.  (13)  is  the  sum  and 
difference  of  three  third  order  polynomials  in  T,  the  e^  and  c^ 
terms  can  be  either  positive  or  negative  depending  of  the  rela- 
tive values  of  m,  n,  r ,  q,  and  s.   Therefore,  no  definite  conclu- 
sion as  to  the  number  of  zeros  and  relative  extrema  of  Eq .  (14) 
can  be  made.   Several  practical  examples  were  solved  and  in  all 
cases  the  transfer  function  had  only  one  extremum.   Since  no 
practical  network  constants  were  found  which  would  yield  a  more 
general  transfer  function  this  network  was  abandoned. 

The  second  basic  network  examined  was  of  the  Wheatstone 
bridge  form.   The  general  form  of  this  network  is  given  in  Fig. 
10.   The  voltage  transfer  function  of  this  network  using  the 
linear  approximation  is 
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Fig.  10a   A  single  section  bridge  network  with 
one  active  element. 


nT 


Fig.  10b   A  single  section  bridge  network  with 
two  active  elements. 
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V 

o    R  -  m  +  nT  (15) 

V   =  2R  +  2m  -  2nT 

At  the  temperature  of  the  bridge  balance  point  the  transfer  ratio 
has  the  value  zero.   The  balance  point  may  be  controlled  entirely 
by  the  values  of  the  fixed  bridge  arms.   Hence,  the  location  of 
the  transfer  function  zero  may  be  determined  by  the  values  of  the 
fixed  bridge  arms.   A  discussion  of  bridge  theory  by  Stout  (2) 
indicates  that  the  sensitivity  of  the  bridge  to  temperature  changes 
can  be  increased  by  placing  two  identical  thermistor  elements 
in  opposite  bridge  arms  as  illustrated  by  Fig.  10b.   This  method 
of  increasing  sensitivity  is  employed  in  the  remainder  of  this 
report.   Two  arrangements  of  the  thermistor  elements  in  the 
bridge  arms  are  possible.   They  give  the  two  transfer  functions 
plotted  in  Fig.  11a  and  lib.   The  curve  of  Fig.  11a  is  the 
negative  of  the  curve  in  Fig.  lib.   The  transfer  function  is  not 
linear,  but  it  is  nearly  enough  to  allow  the  following  intuitive 
reasoning.   Two  sections  with  different  zero  locations  may  be 
cascaded.   Neglecting  loading  effects,  the  combined  transfer 
function  should  be  a  function  in  T  having  two  zeros  as  shown  in 
Fig.  12.   If  three  sections  are  cascaded,  the  output  will  have 
three  zeros  in  T  with  two  relative  extrema.   The  three  zeros  may 
be  located  anywhere  in  the  temperature  range  by  changing  the 
balance  points  of  the  individual  sections.   By  using  the  shifting 
technique,  a  curve  may  be  produced  which  has  one,  two,  or  three 
zeros  in  the  range  of  the  linear  approximation.   It  seems  that 
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Temperature 


Fig.  11a 


Temperature 

Fig.  12b   Two  possible  transfer  functions  for 
the  single  section  bridge. 
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as  many  sections  as  desired  could  be  cascaded,  giving  a  curve 
having  any  given  number  of  zeros  in  T. 

Voltage  Transfer  Function  for  the  Bridge  Network 


The  next  step  was  to  verify  this  intuitive  reasoning  by 
computing  the  actual  transfer  function  for  a  network  consisting 
of  three  cascaded  bridge  sections.   The  transfer  function  should 
include  loading  effects  and  use  practical  values  of  thermistor 
constants  and  bridge  arm  values.   The  circuit  diagram  of  a  three 
section  bridge  is  given  in  Fig.  13.   Due  to  the  complexity  of  the 
node  and  loop  equations  for  this  network,  the  topological  formulas 
of  Seshu  and  Reed  (3)  were  used  to  compute  the  transfer  function. 
These  formulas  relate  the  transfer  function  to  the  structure  of 
the  circuit.   A  major  advantage  of  this  method  is  that  there  is 
no  cancellation  of  terms  as  in  node  or  loop  analysis. 

In  order  to  give  the  topological  formulas,  several  defini- 
tions pertaining  to  the  network  graph  must  be  given.   A  tree  of 
a  graph  is  defined  as  a  connected  subgraph  containing  all  the 
vertices  of  the  original  graph  but  no  closed  circuits.   A  2-tree 
is  a  pair  of  unconnected  subgraphs,  each  subgraph  itself  being 
connected,  containing  all  the  vertices  of  the  original  graph  but 
having  one  less  edge  than  a  tree.   A  2-tree  is  represented  by  the 
symbol  (2-T) .   A  2-tree  admittance  product  is  the  product  formed 
by  multiplying  the  admittances  represented  by  the  edges  of  a 
2-tree.   A  2-tree  admittance  product  is  represented  by  the  symbol 
W.   Consider  the  generalized  graph  represented  by  Fig.  14.   A 
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Fig.    14      Generalized    two-port   network    graph, 
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number  of  2-trees  of  this  graph  will  have  vertices  1  and  2  in 
separate  subgraphs.   Such  2-trees  are  represented  by  the  sym- 
bology  (2-T)   ,.   Other  2-trees  will  have  1  and  2  in  separate 
subgraphs  and  will  also  have  3  and  4  in  separate  subgraphs. 
These  are  represented  by  (2-T)13>24  and  (2-T)14,23'  where  (2_T)13,24 
denotes  those  2-trees  such  that  vertices  1  and  3  are  in  one 
connected  subgraph  and  2  and  4  are  in  the  other  subgraph  and 
(2-T)  ,  ,-  denotes  the  remaining  2-trees.   The  corresponding 
2-tree  admittance  products  are  represented  by  W1>2>  w13i24'  and 

W      .   The  sum  of  all  2-tree  admittance  products  of  a  specified 
14  ,23 

class  is  represented  by  EW.   The  topological  formula  for  the 
voltage  transfer  function  of  a  two  port  is 

v     iw,,  ,,  -  rw  .  ,, 

out  =    13,24 14.23  (16) 

V.  ZW,     9 

in  1 ,2 

Formulas  for  the  number  of  2-trees  of  a  particular  graph  and  the 
method  of  finding  those  2-trees  are  rather  complicated.   The 
problem  can  be  reduced  to  the  much  easier  problem  of  finding 
the  trees  of  a  simillar  graph.   If  terminals  1  and  2  of  Fig.  14 
are  shorted  together  and  the  trees  of  this  graph  are  found,  these 
trees  will  be  the  same  as  the  2-trees  of  the  original  graph  having 
1  and  2  in  separate  subgraphs.   Likewise  if  3  and  4  are  shorted, 
the  2-trees  of  the  original  graph  having  3  and  4  in  separate  sub- 
graphs can  be  found.   Those  2-trees  represented  by  (2-T)13>24  and 
(2-T)    ,   will  be  those  2-trees  common  to  both  the  (2-T)^,  2  and 
the  (2-T)3  4  lists. 


31 


Since  the  method  of  determining  the  trees  of  a  large  graph 
is  rather  involved,  it  is  very  helpful  to  know  the  actual  number 
of  trees  before  attempting  to  determine  them.   For  a  large  graph 
this  is  necessary  as  a  check  to  make  certain  that  all  the  trees 
have  been  found.   The  number  of  trees  of  a  graph  may  be  found  by 
evaluating  the  determinant  of  the  following  matrix, 


[E]  ■  [e.j] 


(17) 


.th 


where  row  i  and  column  i  corresponds  to  the  i    vertex  of  the 
graph,  i  =  1,  2,  3,---,n-l  (for  a  graph  of  n  vertices).   The 
element  e    is  the  number  of  edges  incident  at  the  i    vertex. 
The  element  -e    =  -e..,  i  1    j  ,  is  the  number  of  edges  connecting 
the  1    and  j    verticies.   The  E  matrix  for  the  three  section 
bridge  of  Fig.  14  with  the  input  shorted  is 


[E]  = 


4 

-2 

0 

0 

0 

0 

2 

4 

-1 

-1 

0 

0 

0 

-1 

4 

0 

-1 

-1 

0 

-1 

0 

4 

-1 

-1 

0 

0 

-1 

-1 

2 

0 

0 

0 

-1 

-1 

0 

2 

(18) 


The  determinant  of  E  is  equal  to  256.   Since  the  circuit  is 
symmetrical  there  are  256  (2-T)1  2  and  256  (2-T)3  ^.   The  trees 
of  the  shorted  graphs  may  be  found  by  a  method  outlined  by 
Percival  (4)  and  Illustrated  by  the  example  below.   Using  as  an 
example  the  graph  of  Fig.  15a,  select  two  terminals,  for  example 


32 


Fig.  15a 


Fig.  15b 


Fig.  15c 


Example  of  finding  the  trees  of  a  graph  by 
Percival's  method. 
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2  and  4,  and  construct  the  two  graphs  as  shown  in  Fig.  15b  and 
15c.   The  trees  of  Fig.  15b  are  all  of  the  trees  of  Fig.  15a 
which  contain  the  edge  c;  and  the  trees  of  Fig.  15c  are  all  of 
the  trees  of  Fig.  15a  which  do  not  contain  the  edge  c.   This 
process  is  continued  until  the  trees  of  the  individual  graphs  so 
created  may  be  written  by  inspection.   The  union  of  the  sets  of 
trees  of  these  individual  graphs  is  equal  to  the  set  of  all 
trees  of  the  original  graph.   The  trees  of  Fig.  15a  are  as  fol- 
lows :  c(a+b+d) (e+f )+abde+abdf+adef+abef+bdef =ace+acf+bce+bcf+ 
dce+dcf+abdf+abde+adef+abef+bdef. 

The  voltage  transfer  function  of  the  network  of  Fig.  13  was 
found  by  compiling  the  lists  of  all  appropriate  2-trees  and  re- 
placing each  edge  by  the  admittance  of  that  arm  to  obtain  the  2- 
tree  admittance  products.   Then  the  proper  sums  of  2-tree  admittance 
products  were  obtained  to  satisfy  Eq.  16.   It  seemed  possible  at 
this  point  that  by  using  the  linear  approximation  for  the  ther- 
mistor element  and  leaving  the  fixed  resistor  as  unknowns  a 
direct  expression  in  powers  of  T  could  be  obtained.   Then  relating 
this  expression  to  f(t)  the  fixed  resistor  unknowns  could  be 
determined.   However  the  complexity  of  this  approach  due  to  the 
large  number  of  terms  involved,  256  in  the  denominator,  led  the 
author  to  abandon  this  idea.   Instead  the  computer  program  of 
Appendix  D  was  written,  incorporating  the  program  of  Appendix  B. 
This  is  an  iterative  program  computing  the  values  of  all  2-tree 
admittance  products  for  a  given  temperature  and  summing  them  to 
obtain  the  voltage  transfer  ratio  for  that  temperature.   The 
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value  of  T  is  increased  and  the  process  is  then  repeated  until 
the  transfer  function  for  the  desired  temperature  range  has  been 
found.  This  program  uses  actual  thermistor  constants  and  bridge 
arm  values.  The  only  factor  not  included  in  the  program  was  the 
possible  effects  of  thermistor  self -heat in g .  Previous  intuitive 
reasoning  seems  to  be  well  substantiated  by  the  results  of  the 
program  for  a  sample  problem  as  compiled  in  Appendix  D,  Table  1. 

The  attentuation  of  the  three  section  bridge  is  rather  high. 
If  in  the  example  of  Appendix  D,  10  volts  were  applied  to  the  in- 
put of  the  bridge  the  maximum  output  voltage  would  be  0.09  volts. 
To  be  of  practical  use  this  voltage  would  need  to  be  amplified. 
This  could  be  done  by  using  a  low  gain,  high  stability  d.c. 
operational  amplifier  or  by  chopping  the  output  and  using  a 
higher  gain  a.c.  amplifier.   The  zero  locations  of  the  output 
are  controlled  by  the  fixed  bridge  arms.   Additional  flexibility 
could  be  achieved  by  using  the  voltage  shifting  principle  dis- 
cussed previously.   Therefore  the  three  section  cascaded  bridge 
network  seems  to  have  sufficient  flexibility  to  be  used  as  a 
general  compensation  network  for  the  generation  of  a  variety  of 
temperature  compensation  voltages. 
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COMPENSATING  A  CRYSTAL  RESONATOR 

The  voltage  generation  method  of  the  last  chapter  was  tested 
by  attempting  to  generate  a  control  voltage  for  an  actual  tempera- 
ture compensation  problem.   The  problem  chosen  was  the  elimination 
of  temperature  induced  frequency  variations  from  a  quartz  crystal 
resonator. 

The  frequency  of  a  quartz  crystal  is  rather  sensitive  to 
variations  in  temperature.   This  variation  was  expressed  by 
Bechmann  (5)  as  a  third  order  function  of  temperature.   The  fre- 
quency deviation  is 

Af   _   _   /l.   T.   1   j_  _   /T_T   1  ^   _,_  _   fT_T   1  ^  (19) 


!f.    mjd-V  +  m2(T-To)'  +  m3(T-To)J 


where  m. ,  m, ,  and  m,  are  coefficients  which  depend  upon  the  angle 
of  cut  of  the  crystal,  T  is  the  temperature  of  interest,  Tq  is 
some  reference  temperature,  f  is  the  crystal  frequency  at  tempera- 
ture T  ,  and  if  is  the  deviation  of  crystal  frequency  from  f. 
o 

This  frequency  deviation  may  range  from  as  little  as  +  10  ppm  to 
+  70  ppm  or  more  over  the  temperature  range  -40°  C  to  +70"  C. 
Modern  communication  procedure  requires  a  tolerance  of  not  more 
than  +  1  or  2  ppm  on  resonator  frequency.   In  stationary  operations 
this  deviation  may  be  eliminated  by  placing  the  crystal  in  a 
thermostatically  controlled,  electrically  heated  oven.   The  cry- 
stal temperature  remains  constant  and  therefore  the  deviation  is 
eliminated.   The  power  required  by  the  oven  and  the  weight  of 
oven  insulation  make  this  method  impractical  for  mobile  operations. 
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A  method  of  compensating 

a  crystal  for  mobile  operations  was 

described  by  Malik  (6). 

This  method  will  be  used  with  basic 

alterations.   First,  the 

three  section  bridge  will  be  used  for  the 

generation  of  the  compensating  voltage.   This  is  a  more  complex 

and  versatile  network  an 

d  should  produce  a  better  compensation. 

Second,  computer  techniq 

ues  were  used  allowing  the  constants  of 

the  compensation  network 

to  be  determined  more  easily. 

The  crystal  frequency  may  be  varied  by  changing  the  load 

capacitance  of  the  cryst 

al.   The  variation  is 

y  _   M  +  c° 

x  106 

CO) 

f  -  "M  +  2r 

(C  +C  )     '                          (      ' 

0    X 

where  M  is  a  constant  de 

termined  by  the  value  of  C   desired  at 

some  frequency  f,  C   is 

the  parallel  capacitance  of  the  crystal, 

r  is  the  ratio  of  series 

to  parallel  crystal  capacitance,  and  C^ 

is  the  load  capacitance. 

If  the  capacitance  C   is  varied  by  the 

f                X 

proper  amount  at  each  te 

mperature,  then  the  frequency  deviation 

will  be  eliminated.   The 

variable  capacitance  can  be  provided  by 

a  varactor  diode  placed 

across  the  crystal  and  biased  by  the 

proper  d.c.  control  volt 

age.   The  compensation  network  is  shown 

in  Fig.  16.   The  two  varactorsare  placed  in  series  to  avoid  the 

possibility  of  their  being  self-biased  by  the  a.c.  signal.   The 

parallel  capacitance  Cf 

enables  the  varactor  to  be  operated  in 

the  most  desirable  part 

of  its  characteristic  curve.   The  varactor 

diode  capacitance  is 

C 

-  \            .                              <21> 

vA 

.-1 

m       « 

It           u 
U             01 

o         >i 
3           W 

37 

, 

. 

u 

A/ 

\ 

AAA 

•H 

3 

V   V 

'  V  V 

c£ 

O 
U 

O 

c 

0 

•H 

2 

\                H 

S 

01 

0 

01 

■u 

^            01 

c 

o 

J »        to 

oi 

ft 

R           B 

c 

■o 

0          *J 

B 

O         rH 

o 

IM*-* 

(-J 

■                     n                             t 

u 

(J 

o 

u          > 

w 

o 

•o 

B 

eg 

W 

P 

01 

BJ 

>. 

> 

M 

3 

[ 

M 

■H 

fc 

A    / 

\ 

AA/i 

v  V 

'  V  V 

eJ 

4J 

•H 

3 

O 

U 

•H 

M  ' 

O 

kl 

o 

4J 

> 

<.                                 B 

r-l                                                     V 

W 

•H 

U 

oi 

O 

38 


where  K  and  A  are  constants  and  V  is  the  bias  voltage.   The  net 
capacitance  of  the  two  series  varactors  is 


(22) 


while  the  total  load  capacitance  is 


C   =  C,  +  C 
x     r     £ 


(23) 


The  solution  of  Eq.  (20),  (21),  (22),  and  (23)  give  the  necessary 
d.c.  voltage  needed  to  correct  a  particular  frequency  deviation. 
Starting  with  a  plot  of  frequency  deviation  vs.  temperature  for  a 
crystal,  the  control  voltage  necessary  to  compensate  the  crystal 
at  each  temperature  can  be  computed.   This  is  the  voltage  which 
must  be  supplied  by  the  output  of  the  three  section  bridge, 
suitably  amplified  and  shifted.   The  temperature-voltage  relation- 
ship necessary  for  this  compensation  method  may  be  found  by  three 
different  methods.   The  voltage  may  be  found  by  graphically  pro- 
jecting each  point  around  the  curves  of  Fig.  17.   The  same  data 
may  be  obtained  by  the  solution  of  the  equations  by  the  computer 
program  of  Appendix  E.   Both  of  these  methods  allow  the  projection 
of  frequency  tolerance  curves  into  voltage  tolerance  curves.   A 
more  accurate  method  is  to  actually  build  the  resonator  circuit 
with  the  varactor  diodes  and  fixed  capacitor.   Then  this  whole 
circuit  is  placed  in  a  temperature  variable  oven.   The  frequency 
of  the  resonator  is  measured  by  a  frequency  meter  while  the  con- 
trol voltage  is  supplied  by  a  variable  d.c.  supply  voltage.   As 
the  temperature  is  varied  a  plot  of  frequency  deviation  and 
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control  voltage  necessary  to  correct  the  frequency  is  obtained. 
This  method  is  best,  since  it  relies,  not  on  theory,  but  on 
actual  circuit  operation.   Therefore  it  is  generally  used  for 
practical  problem  solutions. 

Since  a  voltage  curve  obtained  by  the  experimental  method 
was  not  available,  the  frequency-temperature  data  of  Table  1, 
Appendix  E  was  used  in  conjunction  with  the  computer  program  of 
Appendix  E  to  obtain  a  voltage- temperature  curve.   This  data  is 
compiled  in  Table  2,  Appendix  E.   The  frequency- temperature  data 
of  Table  1  is  actual  data  for  a  Midland  one  megacycle  crystal, 
however  crystal  constants  were  not  available.   Therefore,  reason- 
able crystal  constants  were  assumed.   The  voltage-temperature 
relationship  to  be  obtained  as  the  output  of  the  control  network 
is  given  in  Fig.  18a.   The  temperature  range  of  this  plot  is  -40 
to  +70°  C.   This  is  a  110°  C  range,  very  close  to  the  total 
linearized  range  of  the  thermistor.   The  output  of  the  program 
for  the  three  section  bridge  is  in  terms  of  a  voltage  transfer 
ratio,  while  the  information  desired  is  in  the  form  of  voltage 
data.   The  problem  can  be  solved  simply  by  considering  the  data 
from  the  three  section  bridge  program  as  the  voltage  output  of  a 
bridge  energized  by  a  1  volt  input. 

Since  the  bridge  output  voltage  may  be  shifted  at  will,  the 
desired  voltage  curve,  Fig.  18a  may  be  likewise  shifted  as  an  aid 
to  synthesis.   The  curve  of  Fig.  18a  has  two  relative  extrema  and 
no  zeros.   This  curve  was  shifted  to  have  a  zero  midway  between 
the  two  extrema  as  shown  in  Fig.  18b.   As  an  aid  to  the  location 
of  the  other  two  zeros  the  piecewise  linear  approximation  and  the 
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Temperature  ° C 
Fig.  18a   Desired  control  circuit  output 
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18b   Desired  control  circuit  output  shifted 
to  provide  a  zero  location. 
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idea  of  cascading  without  loading  effects  were  used.   The  output 
of  a  single  bridge  section  containing  as  its  temperature  variable 
element  the  piecewise  linear  approximation  is  given  in  Fig.  19. 
If  three  of  these  sections  were  cascaded  to  obtain  the  output  of 
Fig.  18b,  the  zero  locations  would  need  to  be  as  shown  in  Fig.  20. 
These  zero  locations  are  192°  K,  301°  K,  and  438°  K  and  are  used 
as  the  first  trial  in  the  solution  of  the  transfer  function  pro- 
gram. 

In  order  to  implement  the  program,  the  fixed  bridge  arm 
values  corresponding  to  the  zero  locations  above  must  be  found. 
After  picking  the  desired  thermistor,  the  values  of  the  parallel 
thermistor-resistor  combination  are  computed  for  values  from  150 
K  to  450°  K  and  compiled  as  Table  3,  Appendix  E.   To  have  a  zero 
at  192°  K  the  fixed  resistor  of  the  first  bridge  section  must  be 
equal  to  the  parallel  thermistor-resistor  resistance  or  9,964 
ohms.   Using  the  data  from  Table  3,  Appendix  E  the  fixed  resistor 
values  for  the  three  section  were  found  to  be  9,964,  4920,  and 
410  ohms.   The  thermistor  chosen  had  constants  of  B  =  3,000°  K, 

I   =  300°  K.  and  R   =  10,000  ohms.   This  data  was  entered  into  the 
o  '       o 

bridge  transfer  function  program.   The  first  trial  did  not  yield 
satisfactory  results.   The  zero  locations  were  amended  as  seemed 
necessary  and  several  more  trials  made.   The  first  concern  was  to 
locate  the  relative  extremum  at  the  proper  temperature  then  to 
adjust  the  general  shape  of  the  curve.   When  adjusting  the  ex- 
trema  locations  the  higher  and  lower  temperature  zeros  were  moved 
independently.   Moving  a  zero  location  by  10  degrees  would  result 
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in  a  5  degree  movement  in  the  location  of  the  adjacent  extremum. 

After  the  extrema  were  located,  adjusting  the  curve  shape  required 

the  simultaneous  adjustment  of  all  three  zero  location  in  order 

to  preserve  the  locations  of  the  relative  extrema.   After  several 

trials  a  curve  of  satisfactory  shape  was  obtained  except  that  it 

was  shifted  toward  the  higher  temperatures  as  shown  in  Fig.  21. 

The  value  of  T   was  changed  to  T   =  273°  K.   The  information  of 

0                          0 

Table  3  was  recomputed  using  the  new  thermistor  data.   Appropri- 

ate zero  location  and  hence  resistor  values  were  chosen.   The 

bridge  output  voltage  curve  of  this  trial  was  located  properly 

on  the  temperature  axis  but  the  curve  was  too  peaked  at  the  ex- 

trema as  illustrated  by  Fig.  22.   Therefore  the  value  of  B  was 

changed  to  B  =  2,500°  K  and  again  the  thermistor-resistor 

parallel  resistance  was  computed.   Values  of  the  fixed  resistors 

of  8,860,  3,290,  and  845  ohms  were  chosen  for  the  three  bridge 

sections  representing  zero  locations  of  233,296,  and  369  degrees 

Kelvin.   The  bridge  transfer  function  was  computed  and  seemed 

to  match  rather  well  the  voltage  curve  of  Fig.  18a.   This  trans- 

fer function  data  is  compiled  in  Table  5,  Appendix  E.   The 

bridge  transfer  voltage  was  then  amplified  and  shifted  to  give 

the  actual  control  voltage. 

The  input  voltage  of  the  three  section  bridge  was  taken  to 

be  one  volt.   This  value  was  chosen  to  satisfy  the  current 

limitations  of  the  thermistor.   The  maximum  current  would  be 

about  0.3  milliamperes  which  should  not  cause  undue  thermistor 

heating.   The  voltage  amplification  factor  was  obtained  by 
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Temperature 

Fig.  21   A  control  circuit  output  function 
showing  the  effects  of  too  high 
a  value  of  T  . 


Temperature 

Fig.  22   A  control  circuit  output  showing 

the  effects  of  too  high  a  value  of 
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Amplification  _  Max  desired  voltage  -  Min  desired  voltage    -25j 
Factor         Max  output  voltage  -  min  output  voltage 

Using  Eq.  (25)  the  necessary  amplification  was  computed  to  be 


Amplification    8.5350  -  6.7985 
Factor       "  .01233  -  .00725 


88.678 


(26) 


The  necessary  shifting  is 


Shifting  m    Hax  deslre<j  voltage  -  (Amplification  x 


Voltage 


Max.  output  voltage) 


(27) 


It  was  necessary  to  reverse  the  polarity  of  the  bridge  output 
voltage,  however  in  an  actual  circuit  this  would  only  require  the 
reversal  of  the  connection  terminals.   For  this  problem, 


Shifting  Voltage  =  8.5350  -  88.678  x. 01233  -  7.4715 


(28) 


Applying  the  amplification  and  shifting  of  Eq.  (26)  and  (28)  to 
the  data  of  Table  5,  Appendix  E,  the  actual  control  voltage  of 
Table  6,  Appendix  E  was  obtained.  This  information  was  plotted 
in  Fig.  23  along  with  the  desired  control  voltage  curve  and  the 
tolerance  voltages  which  will  allow  a  +  1  ppm  deviation  in  reson- 
ator frequency.  The  control  voltage  is  within  these  tolerances 
from  -35  to  +70°  C. 

A  study  of  allowable  tolerances  on  amplifier  gain  and  on 
shifting  voltage  was  not  attempted.   However  during  the  process 
of  determining  the  proper  amplification  and  shifting  it  was  noted 
that  a  variation  of  up  to  6  percent  in  the  amplification  factor 
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would  cause  only  about  1  ppra  degradation  in  the  frequency 
tolerances.   Similiarly  it  was  noted  that  the  shifting  voltage 
would  need  to  be  held  within  about  1  percent. 

The  network  shown  in  Fig.  24  will  generate  a  control,  vol- 
tage which  will  compensate  a  Midland  one  megacycle  crystal  placed 
in  the  circuit  of  Fig.  17  to  within  +  1  ppm  over  a  105"  C  tempera- 
ture range.   Previously  described  methods  of  voltage  generation 
(6)  had  yielded  a  compensation  of  approximately  +  2.5  ppm. 
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CONCLUSION 

Certain 

temperature  compensation  problems  require  the 

generation  of 

a  temperature  compensation  voltage.   Various  pro- 

blems  and  considerations  of  designing  thermistor  networks  were 

discussed  in 

this  report.   The  cascaded  bridge  network  was  de- 

veloped  and  c 

emonstrated  to  be  a  versatile  method  of  generating 

temperature  compensated  voltages.   This  network  along  with  the 

voltage  shift 

ing  technique  was  applied  to  the  problem  of  compen- 

sating  a  crys 

tal  resonator  for  temperature  induced  frequency 

variations . 

The  results  of  theoretical  computations  using  this 

network  show 

that  frequency  deviation  can  be  reduced  to  +  1  ppm 

or  less  .   Coir 

puter  techniques  as  listed  in  the  Appendicies  were 

developed  to 

aid  in  the  analysis  of  this  network  and  in  the 

determination 

of  various  network  constants. 

Several 

methods  of  increasing  the  versatility  of  the  cas- 

caded  bridge 

circuit  became  evident  in  the  latter  stages  of  this 

investigation 

but  time  was  not  available  to  examine  them  fully. 

The  performance  of  one  stage  of  the  bridge  is  affected  by  the  re- 

sistive load 

applied  to  its  input  and  output  terminals  by  the  pre- 

ceeding  and  f 

ollowing  stages.   This  effect  may  be  varied  by 

altering  the 

impedance  level  of  the  various  stages.   Two  transfer 

functions  were  computed;  in  one  the  impedance  level  was  increased 

by  a  factor  c 

f  ten  per  section,  in  the  other  the  impedance  level 

of  all  stages 

was  the  same.   The  effects  of  this  impedance  change 

were  not  great,  however  a  proper  understanding  of  loading  effects 

might  enable 

the  designer  to  use  them  in  a  limited  but  beneficial 

manner . 
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In  this  report  the  same  thermistor  was  used  in  each  section 
of  the  bridge,  only  the  fixed  arms  were  changed.   This  aids  the 
simplicity  of  the  design  and  would  contribute  to  manufacturing 
economy.   As  stated  previously  the  change  of  the  B  constant  of 
the  thermistor  affected  the  slope  of  the  linearized  section.   The 
change  of  the  B  or  T   constants  of  the  thermistor  in  only  one 
section  of  the  bridge  might  enable  the  designer  to  alter  the  shape 
of  only  a  portion  of  the  transfer  curve. 

The  zero  locations  of  the  transfer  function  are  determined 
by  the  fixed  arm  values  of  the  various  bridge  sections.   In  this 
study  the  first  or  low  temperature  zero  location  was  provided  by 
the  first  bridge  section,  the  second  zero  by  the  second  bridge 
section,  and  the  high  temperature  zero  by  the  third  bridge  sec- 
tion.  There  is  no  reason  that  this  order  could  not  be  reversed 
or  scrambled.   For  example,  the  high  temperature  zero  could  be 
provided  as  a  result  of  the  first  section,  the  low  temperature 
zero  by  the  second  section,  and  so  forth.   Experience  gained  as  a 
result  of  this  study  led  the  author  to  believe  this  technique 
might  have  significant  effects  on  the  transfer  function  of  the 
bridge.   Although  no  positive  proof  was  available. this  as  well  as 
the  other  factors  mentioned  would  be  worthy  of  further  investigation. 
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Program  1 


C 
3000, 


PR 
B=T 

0  = 
RC  = 
TC  = 
REA 
F  =  E 
G 
G  =  R 
DG 
DO 

R  =  D 

D!  = 

D2 

03 

DA 

D4 

DB 

D5  = 

DC  = 

D6  = 

AO, 

A0= 

Al  = 

A2  = 

A3 

A4  = 

A  5 

A6 

CI 

CI 

C2  = 

C3 

C4= 

C5= 

C6  = 

C7  = 

T  =  2 

■R  = 

R4  = 

R  =  G 

R4 

PUN 

T  =  T 

IF( 

CO" 

END 

TYP 


ISTCR  BY  A  POWER  SERIFS 

STCR 

HE  POWER  SERIES  IS  EXPANDED 

R 

STCR 


GRAM  FOR  APPROXIMATING  A  THFRV 
HF  BETA  CONSTANT  CF  THE  THERMI 
THE  TEMPERATURE  AROUND  WHICH  T 

THE  RO  VALUE  OF  THE  THERMI STC 

THE  TO  CONSTANT  OF  THE  THERMI 
D»B>QiRC»TC 
XPI-B/TC) 

THE  COMBINATION  CF  ALL  CONSTANT  TERMS  CF  THE  THERMISTOR  EQUATION 
0*F 

G*EXP(B/Q) 

D1»...D6=THE  THERMISTOR  EQUATION  AND  ITS  DERIVATIVES  EVALUATED 
AT  THE  0  POINT 


B*R/(Q*Q 
2.*B*R/(Q 
-6.*B*R/( 

2A.*B*R/( 
DA+e**4*R 
120.*B*R/ 

-de-; j.*b 

72C.*B*R/ 

DC+3C0.*B 

A  i  ,  .  .  .  A  6  = 

DO 

Dl 

D2/2. 

D3/6. 

D4/24. 

05/120. 

D6/720. 

C2»...C7= 

A0-A1*Q+A 

A1-2.*A2* 


) 

*Q*Q)+R*B*R/ |Q*Q*Q*Q) 

G*Q*Q*Q) -6.*B*B*R/ (0* 

Q**5)+36.*B*B*R/(Q**6 

/(Q**8 ) 

(U'**6  1+240. *B*B*R/  (  U* 

**4*R/ ( Q**9)-B**5*R/ ( 

IQ**7)+1800.*B*B*R/ (Q 

**4*R/ ( o**i c |+30.*B»* 

THF  RESPECTIVE  D  TERM 


*5)-tS*6*B*R/<Q**6) 
>+12.*B**3*R/(Q**7) 

»7)+120.*o**3*R/(Q**8) 

Q**10) 

**8>+1200.*B*B*d*R/(Q**9 ) 
5*R/(Q**11 )+b**6*R/ (0**12 ) 

S  DIVIDED  BY  THE  PROPER  FACTORIAL 


A2-3.*A3* 
A3-4.*A4* 

A4-5.*A5* 
A5-6.*A6* 

A6 
33. 

THE  ACTUA 
THE  POWE 
*EXP( B/T) 
C1+C2*T+C 
CH.T.R.R4 
+  2. 

T-334.)  3 
TINUE 


THE  POWERS  OF  T  IN  It 

2*a*Q-A3*Q*0*Q+A4*Q*C 

G  +  3  .  *A3*Q*Q-4  .  *A  4*G  *Ci 

G+6.*A4*Q*Q-10.*A5*Q* 

Q+10.*A5*O*Q-20.*A6*C 

Q+l 5.*A6*0*0 

n 


EXPANDED  POWER  SERIES 
*(.*ti-A5*0**5+A6*Q**6 
*d+5.*A5*U'**4-6.*A6*Q**5 

*0+15.*A6*Q*Q*Q*Q 
*Q*0 


L  THERMISTOR  RES  I  STAN 
R  SERIES  APPROXIMATE 


'T*T+C4*T*T*T  +  Cb*T*T 


CE 

N  TC  TnE  THERMISTOR  RESISTANCE 

*T*T+C6*T**5+C7*T**6 


ICAL  INPUT  DATA 
283.0        100C 


300. 


Table    1 
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Temperature  °K 

Thermistor 

Power  Serii 

2.33O0E+02 

1.7735E+05 

1.6882E+05 

2.35jGE+02 

1.5894E+05 

1.5293E+05 

2.370OE+02 

1.4271E+C5 

1.3835E+05 

2.39OCE+02 

1.2837E+05 

1.2530E+05 

2.4100E+02 

1.1567E+05 

1.1339E+05 

2.4300E+G2 

]  .0440E+G5 

1.0284E+C5 

2.450OE+02 

9.4395F+G4 

9.3520E+C4 

2.470OE+02 

8.5484E+04 

8.4990E+04 

2.4900E+02 

7.7539E+04 

7.724GE+G4 

2.510CE+02 

7.0441E+C4 

7.0590E+04 

2.5300E+02 

6.4G9GE+G4 

6.416GE+G4 

2.5500E+02 

5.8399E+C4 

6.6buuE+04 

2.57OQE+02 

5.3289E+04 

5.367ot+U4 

2.59uOE+02 

4.8696E+04 

4  •89^uE  +  l)4 

2.6100E+U2 

4.4b6GE+04 

4.471GE+04 

2.631.CE  +  C2  ' 

4.083CE+04 

4.097OE+04 

2.65O0E+02 

3.7462E+04 

3.7390E+04 

2.6700E+02 

3.4417E+04 

3.4450E+C4 

2.69COE+02 

3.1658E+04 

3.2'JlOE  +  C4 

2.7100E+02 

2.9157E+04 

2.94GGE+04 

2. 73006+02 

2.6886E+C4 

2.7200E+04 

2.75COE+0? 

2.4821E+04 

2.3v9GE+G4 

2.77C0E+Q2 

2.2941E+04 

2.3^lGE+U4 

2.790CE+02 

2.1227E+04 

2.  loVOE+04 

2.81uOE+02 

1.9S63E+G4 

2.0j9ut+04 

2.8300E+02 

1.8234E+04 

1.8560E+04 

2.8500E+O2 

1.6927E+04 

1 .70O0E+04 

2.8700E+02 

1.5730E+04 

1.5780E+04 

2.890CE+02 

1.4632E+04 

1  .51uoE+u4 

2.9100E+C2 

1.3624E+04 

1 .4j30E+04 

2.93006+02 

1.2699E+04 

1.3200E+04 

2.950CE+02 

1.1847E+04 

1.2310E+C4 

2.9700E+02 

1.1063E+Q4 

1.1640E+04 

2.99Q0E+02 

1.0340E+04 

1.0810E+04 

3.0100E+02 

9.6732E+03 

1 .0030E+04 

3.0300E+C i 

9.0573E+03 

9.27GGE+03 

3.O5C0E+02 

O.4880E+03 

8 .6 luut+03 

3.070OE+02 

7.9611E+03 

8.29^JE+U3 

3.0900E+02 

7.4732E+03 

7.73O0E+03 

3.1100E+02 

7.02-9E+03 

7.55uot+G3 

3.130CE+C2 

6.6C12E+03 

6.99O0E+03 

3.1500E+02 

6.2115E+03 

6.7300E+03 

3.1700E+02 

5.8492E+03  . 

6.34U0E+G3 

3.1900E+O2 

5.5123E+03 

6.56LuE+G3 

3.2100E+02 

5.19S5E+C3 

6.13GGE+03 

3.230OE+02 

4.9063E+03 

6.35^0E+03 

3.2500E+02 

4.6337E+C3 

5.840GE+03 

3.2700E+02 

4.3793E+03 

6.29UOE+G3 

3.29UOE+02 

4.1418E+03 

6.13OOE+03 

3.3100E+02 

3.9198E+03 

6.32GOE+03 

3.33GOE+G2 

3.7121E+03 

6 .2400L+U3 
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APPENDIX  B 
Program  2 

C   C   PROGRAM  FCR  THE  RESISTANCE  CF  A  THERMI STCR-kES I STOR  IN  PARALLEL 
C      THERMISTOR  CONSTANTSARE  B.  TO.  AND  RC 

C      RP  IS  THE  PARALLEL  RESISTOR,  RS  IS  THE  SERIES  RESISTOR 
3  READ.B.TO.RC 

READ.RP.RS 
C      A,C,  ANf)  F  ARE  CONSTANTS  RELATIVE  TO  THE  THERMISTOR  VALUE 

OB/TO 

T*151. 

i  a=e/t 

D  =  A-C 

F=EXI  (D) 
C      R  =  THE  THERMISTOR  RESISTANCE  AND  RN=  THE  COMBINATION  RESISTANCE 

R=RC*F 

RN  =  (RKRR)'/  (R  +  RPJ+RS 

PUNCH, T.RMtR 

T=T+2. 

[FIT-451.)  1,2,2 
2  CCHTINUE 

GO  TO  3 

END 
C      TYPICAL  INFJT  DATA 

300O;c     3CO.C    1000C.0 
lOOOC.C    O.o 

2500.0     273. 'j      10000.0 
10000.C    O.o 
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APPENDIX  C 
Program  3 


C 


C 

22 
9.0 
XxC 
2 
2 
2 
2 
2 
2 
2 
2 
2 
3 


PLO 

DIM 

]  FCR 

?  FCR 

3  FCR 

4  FCR 
9  RFA 
:  RFA 

SCA 
REA 

6  DC 
15  P(  I 

PUN 
N  =  2 
L  =  l 
DC 
REA 
DC 

7  P(  I 
L»C 
N  =  l 
DC 
M=( 
IF( 
L=M 
IF( 
N=M 

8  IF( 

12  PI!'1 
GC 
P(M 
CON 

13  PUN 
G 
END 
TYP 


TTING  ROUTINE  FCR  FCRGC 

ENSICM  P( 1^8 ) .Al 36) ,B(36 I 

MATIE1C.3.2X.10BA1  ) 

MATI36A2) 

MATI  12X.108A1 ) 

MAT (Z 12) 

D    h  >  <  ,  J 

D    ,  A  MAX, A MTN 

LE=ABS( ( AM AX -AM  IN ) /107. 0 ) 

D    2,     (B< I ) » I=1>K) 

15    1=1.108 

)=B(1) 

Ch    3i(P( I ) » 1=1,106) 

08 

13  NJK=1»J 

D ,     t  A  (  I  )  ,  I  =  1  ,  K  ) 

7    I  =  \',L 

1=0.0 

OS 

14  1=2, < 
A{ I  l-AMINf/SCALE+2.0 

L-M)  5,16,16 


5 
16 
17 


1] 

14 


N-M )  8,8,17 

P ( M ) )  11,12,11 

)=R(  I  ) 

TC  14 

)=.14E-36 

TINUE 

CH  1 ,A( 1 ) , (P( I ) , I=1,L) 

TC  9 

ICAL  INPUT  DATA 


0    6.000 

.330CE+02 
.AIC-OE+02 
.4900E+02 
.57O0E+02 
.65O0E+02 
.7300E+02 
•810GE+02 
.890GE+02 
.970CE+02 
.050CE+C2 


7.2077 

7.0475 
6.9085 
6.8156 
6.7935 
6.8572 
7.0103 
7.241.4 
7.5236 
7.6^63 
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*-SfMICKELSCN 
» 1 0  G  6 


D] 

FC 

F 

FT 

F 

M 

Al 

bC 

RP 

RE 

RE 

A 

A 

D 

F  = 

G  = 

R  = 

Y  = 

CC 

DC 

AM 

AA 

DD 

FF 

EE 

Ul 

B, 

Ul 

U2 

U3 

VI 

B  = 

C  = 

V2 

E* 

H 

V3 

S 

T 

WA 

WA 

WB 

W= 

TH 

Z< 

Z( 


MENS 
RMAT 

RMAT 
RMAT 
RMAT 
»N?> 
.Fli 
.  TO 

AND 
AD  9 
AD  8 
D»F 
l./A 
A 

l./F 
F 

]  ./R 
R 
=  60/ 

2  M 

M 
=  50/ 

AA- 

EXP 
=  RC* 

(EE 
C»E 
=  TH 
=  U1 
=  U1 
=  1./ 
VI 
VI 

1 

V2 
V2 
1./ 

V3 
V3 

AMD 
*A»D 
=  A*D 
WA-W 

Z 
001) 
002) 


ION  Z( 

(F10. 
(2F10 
1 5  F 1C 

(313) 
AND  N3 

AND  R 
.  RC  A 

RS  AR 

N1.N2 
.  A 1  ,  F  1 
G»R.  A 
1 

1 

1 

TO 


APPENDIX  D 
Program  4 


256) 

5) 

.5) 

.1  ) 

DEFINE  THE  TEMPERATURE  RANGE 
1  ARE  THE  FIXED  RESISTORS  OF  THE  BRIDGE  SECTIONS 
RE  THERMISTOR  CONSTANTS 
E  THE  LINEARIZING  RESISTORS 
.N3 

,R1»BC»TC»RCjRP»RS 
ND  Y  ARE  ADMITTANCES  OF  FIXED  ARMS 


=N1.N2.N3 

AM 

CC 

(DD) 

FF 

*RP)/( EE+RPl+RS 

H,S.  AND  T  ARE  THE  ADM1TT; 

E  LINEARIZED  THERMISTOR 


Ul 


NCES  OF  THE  THERMISTOR  ARMS 


/U2 


U3 


WB  ARE  THE  TERMS  OF  THE  TRANSFER  FUNCTION  NUMERATOR 
«F*G«S*T  +  B*C*E*H*S*T  +  A*-D*Irj-'H*R!!-Y  +  B-s-C*-F*G*RR-Y 
*F»G*R*Y+B*C*F*G*S*T+A*D#E*H*S*T+e*C*E*H*R*Y 
B 

TERMS  ARE  THE  DENOMINATOR  OF  TnE  TRANSFER  FUNCTION 
=A*D*F*G*R*Y 
=  A*D*F*G--  R*S 
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Z  (Uu3 

=A#D*F*G*T*Y 

Z1C1.4 

=3 >F*G*R*S*T 

Z  ( 00  5 

=A*D*F*G*S*T 

ZIC06 

=A*D*G*R*S*Y 

ZIC07 

=A*D*F*R*S*Y 

Z(C08 

=A*D*G*R*T*Y 

Z<009 

=A*D*F*R*T*Y 

Z(C1C 

=A*D*G*R*S*T 

Z  ( C 1 1 

=A*D*F*R*S*T 

Z « C 1 2 

=A*D*G*H*R*Y 

Z  t  C 1 3 

=A*D*G*H*R*S 

Z  (  0 1 4 

=A*D*G*HJ»T*Y 

Z  { 0 1 5 

=A*D*G*S*T*Y 

Z  (016 

=A*D*F*S*T*Y 

Z  <  0 1 7 

=A*D-*G*H*S*T 

Z  (  0  1  6 

=A*DfE*F*R*Y 

ZI119 

«A*D*E*F*R*S 

Z  ( 02  0 

=A*D*E*F#T*Y 

Z(021 

=A*D*E*F*S*T 

Z<C22 

=A*D*E*R*T*Y 

Z  (023 

=A*D*E»R*S*Y 

Z  ((.24 

=A*D*H*R#S*Y 

Z  (i-25 

=A*D*H*R*T*Y 

Z(026 

=A*D*E*R*S*T 

Z(C27 

=A*D*H*R*S*T 

Z<^28 

=A*D*E*H*R*Y 

Z(029 

=A*D*E*H*R*S 

Z  (030 

=A*D»E*5*T*Y 

Z  ( 03  ] 

=A*D*H*S*T*Y 

ZJ032 

=A*D*E*H*S*T 

Zt033 

=A*F*G*R*S*Y 

ZI034 

=D*F*G*R*S*Y 

ZIC35 

=A#F*G*R*T*Y 

Z(036 

=A*F*G*S*T*Y 

Z(C37 

=D*F*G*R*T*Y 

Z(C38 

=A*F*G*R*S*T 

ZU'39 

=A*F*G*ri*R*Y 

ZU4G 

=D*F*G*H*R*Y 

ZK;41 

=A*F*G*H*R*S 

ZI042 

=D*F*G*H*R*S 

Z  (C43 

=  A  ;-F*G*H*-T*Y 

Z  ('44 

=D*F*G*H*T*Y 

Z<045 

=D*F*G*S*T*Y 

ZU46 

=A*F*G«H*S*T 

ZIC47 

=A*F*H*R*S*Y 

ZU48 

=0*F*H*R*S*Y 

Z  (049 

=A*F*H*R*T*Y 
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ZC050) 
ZC051! 
Z  ( C  5  2 ) 
Z(053) 
Z(054) 
ZI055) 
ZI056) 
Z  ( C  5  7 ) 
Z (058) 
Z(059) 
21!  60) 
Z(C61 ) 
Z (162) 
Z(C63) 
ZI064) 
Z  CG65) 
Z (C66) 
Z (067) 
ZCC68) 
Z  (<  69) 
ZI07CI 
Z  ( 0  7 1 1 
ZIC72) 
Z(l73) 
Z<074) 
Z(o75> 
Z (076) 
Z  (  0  7  7  ) 
2(078) 
ZI079) 
ZI080) 
Z(081) 
Z  (  0  8  2  ) 
Z  (  0  8  3  ) 
Z(u84> 
2(085) 
Z(086) 
ZI087I 
Z(088) 
ZI089) 
Z (090) 
Z (091 ) 
Z (C97) 
Z (093) 
Z(C94) 
ZIC95) 
ZU'96) 
ZI097) 
Z (098) 


=D*F*H 
=  A*F*H 
=  D*F*H 
=A*F*H 
=D*F*H 
=A  tB*F 
=  A»F»r 
=D*E*F 
=  A«B*F 
=A*E*F 
=A*B*F 
=A*E*F 
=D*E*F 
=d*F*G 
=A*B*F 
=A»E*F 
=D*E*F 
=A*B*G 
=A*E*G 
=D*E*G 

:A*8-"-F 

=A*B*G 
=  A*E*G 
=D*E*G 
=A*d*F 
=A*B*G 
=A*E*G 
=D*E*G 
=  A*B*F 
=A*B*G 
:A*E*G 
=  A*E*F 
=D*E*G 
=D*E*F 
=A*B*G 
=A*E*G 
=A*E*F 
=D*E*F 
=A*B*G 
=  A*E*G 
=D*E*G 
=  A#E«F 
=  D*E*F 
=  A*B*G 
=A  >E*G 
=  D*E*G 
=  A*B*F 
=  A*B*G 
=  A*E*G 


*R*T*Y 
*R*S*T 
*R*S*T 
*S*T*Y 
*S*T*Y 
»G*R*Y 
*G*R*Y 
#G#R*Y 
»G*R*S 
*G*R*S 
*G*T*Y 
*G*T*Y 
»G*T*Y 
*H*S*T 
»G*S*T 
*G*S*T 
*G*S*T 
*R*S*Y 

*  R  *  S  *  Y 
#R*S*Y 
*R*S*Y 
*R*T*Y 
*R*T*Y 
*R*T*Y 
*R#T*Y 
*R*S*T 
*R*S*T 
*R*S*T 
*R*S*T 
*H*R*Y 

H*R*Y 
*H*R*Y 
#H*R*Y 
*H*R*Y 
#h*R*S 
*H*R*S 
*h*R*S 
»H*R*S 
*rt*T*Y 
*H*T*Y 
*H*T*Y 
*H#T*Y 

H*T*Y 
*S*T*Y 

*  s  *  t  *  y 
*S#T*Y 
*S*T*Y 
*H»S*T 

#S*T 
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z 

1.99) 

z 

100) 

z 

101) 

z 

102) 

z 

103) 

z 

1041 

z 

105) 

z 

106) 

z 

107) 

z 

108) 

z 

109) 

z 

110) 

z 

111) 

z 

11?) 

z 

113) 

z 

114) 

z 

1  1?) 

z 

116) 

z 

117) 

7 

1181 

z 

119) 

z 

12  0) 

z 

121) 

z 

122) 

z 

1231 

z 

124) 

z 

123) 

z 

126) 

z 

127) 

z 

128) 

z 

129) 

z 

130) 

z 

13]  ) 

z 

132) 

z 

133) 

z 

134) 

z 

135) 

z 

136) 

z 

137) 

z 

133) 

z 

139) 

z 

140] 

z 

141  ) 

z 

142) 

z 

143) 

z 

144) 

z 

145  ) 

z 

146) 

z 

147! 

z 

143) 

z 

149) 

z 

150) 

=  D*E* 
=A*E* 
=D*E* 
=A*B* 
=  A*B* 

=  A    -RX 

=A*B* 

=A*B* 
=A#B* 
=A*E* 

=  D#E* 
=  A*u* 
=  A*B* 
=  A*E* 
=  D*E* 
=  A*B* 
=A*B* 
=  A  *  E  * 
=D*E* 

=  A-SO* 

= A  *  B  * 
=  A*D* 
=  A*B* 
=A*B* 
=  A*B» 
=A*E* 
=D*E* 
=  A*B* 
=C*D* 
=  C*D* 
=C*D* 
=D#F* 
=  C*D* 
=  C*D* 
=C*D* 

=c*o* 

=  C*D« 
=  C*D* 
=C*D* 
=C*D* 
=C*D* 
=  C«D* 
=  C*D* 
=C*D* 
--C  *D* 
=  C*0* 
=C*D* 
=  C*D* 
=  C*D* 
=  C*D* 
=B*C* 
=  C*D* 


G*H*S*T 
F*H#S*T 
F*H*S*T 
E*F*R*Y 
E*F*R*5 
E#F*T»Y 
E*F#S*T 
E*R#S*Y 
H*R*S*Y 
H*R*S*Y 
H*R*S*Y 
E*R*T*Y 
H*R*T*Y 
H*R*T*Y 
H*R*T*Y 
E#R*S*T 
H*R*S*T 
H*R*S*T 
H*R*S*T 
r*H-*R*Y 

f*h*R*s 
E*H*T*Y 
E*H*T*Y 
E*S*T*Y 
H*S*T*Y 
H*S*T*Y 
H»s#T*Y 
E*n*S*T 
F*G#R*Y 
F*G*R*S 
F»6*T*Y 
G*H*S*T 
F*G*S*T 
G*R*S*Y 
F*R*S*Y 
G*R*T*Y 
F*R*T*Y 
G*R*S*T 
F*R*S*T 
G*H#R*Y 
G*ri*R*S 
G*H*T*Y 
G*S*T*Y 
F*S*T*Y 
G*H*S*T 
E*F*R*Y 
E*F*R*S 
E  *  F  *  T  *  Y 
E*F*S*T 
E*R*S*Y 
H*R*S*Y 
H*R»S*Y 
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Z  (151 ) 
ZU52) 

Z  (  1  5  3  ) 
Z (154) 
Z  (  1  5  5  I 
Z ( 156) 
Z (157) 
Z  (  156) 
Z ( 159) 
Z (160) 
?(161) 
Z (162) 
Z (163) 
Z  (  16'.) 
ZI165) 
Z( 166) 
Z(  167) 
7  ( 168) 
Z1169) 
ZC170] 
Z  (  1  7  1  ) 
ZU72) 
Z(173) 
Z(  174) 
Z(175) 
Z ( 176) 
Z  (  1  7  7  ) 
Z1178) 
ZU79) 
Z(18C) 
ZI181 ) 
Z(182) 
Z( 183) 
Z (184) 
Z (185) 
Z(186) 
Z(] 87) 

zneei 

Z ( 189) 
Z (190) 
Z(191) 
Z(19?) 
ZI193) 
Z (194) 
Z ( 195) 
Z( 196) 
Z  (  197) 
Z  (  198) 
7 (199) 
Z (200) 
ZI201) 
Z (2C2) 
Z(2j3) 


=C*D*E*R*T*Y 
=C*D*H*R*T*Y 
=C*D*F*R*S*T 
=C*0*H*R»S*T 
=  C  *D*E*H*R*Y 
=  C*D*E*ti*R*S 
=  C*D*E*-H*T*Y 
=C*D*E*S*T*Y 
=C*D*H*S*T*Y 
=C*D*E*ri*S*T 
=B*F*G*R*S*Y 
=C*F*G*R*S*Y 
=  RtfF*G«R«T*Y 
=C*F*G«R«T*Y 
=  C*F*G*R-»S*T 
=B*F*G*H*R*y 
=C*F*G*I1*R*Y 
=B*F*G*m*R*S 
=C*F*G*H#R*S 
=B»F*G*H«T*Y 
=  G*F*G---ri*T»Y 
=5*F*G*S*T*Y 
=  C*F*G*S':-T*Y 
=C*F*G*H*S!T 
=B*F*H*R*5*Y 
=C*F*H*R*S*Y 
=6*F*H*R*T*Y 
=C*F*H*R*T*Y 
=B*F*H*K*S*T 
=C*F*H*R*S*T 
=B*F*H*5*T*Y 
=C*F*H*S*1*Y 
=  tj*C*-F*G*R*Y 
=3*E*F*G*R*Y 
=C*E*F*G*R*Y 
=B«C*F*G*R*S 
=n*F*F*G*R*S 
=S*E*F*G*R*S 
=C*E*F*G*R*S 
=B*C*F*G*T*Y 
=o*E«F*G*T*Y 
=C#E*F*G*T*Y 
=8*C*F*G*S*T 
=  0  .-F*G*R*S*T 
=B*E*F*G*S*T 
=C*E*F*G*S*T 
=B*C*G*R*S*Y 
=B*E*G*R*S*Y 
=C*F*G*R*S*Y 
=C*B*F*R»S*Y 
=B*C*G*R*T*Y 
=B*E*G*R*T*Y 
=C*E*G*R*T*Y 
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Z(204 
Z(205 

Z  (2u6 
Z1207 
Z(208 

Z  (2U9 
Z  (  2  1  C 
Z  (211 
Z  (21? 
Z(213 
ZI214 
Z  (215 
Z  (2)6 
ZI217 
Z  (218 
Z  (219 
Z  ( 2  2  0 
Z(221 
Z(222 
Z(223 
Z  (224 
ZI225 
Z(226 
Z  (227 
Z(228 
Z(229 
Z(230 
Z(231 
Z  (232 
Z  (  2  3  3 
Z  (234 
Z  (  2  3  5 
Z(236 
Z(237 
Z(238 
Z  (239 
Z  (240 
Z  (241 
Z(242 
Z(243 
Z  (244 
Z  (245 
Z(246 
Z  (247 
Z  (248 
Z  (249 
Z  ( 2  5  0 
Z(251 
Z(252 
Z  (253 
Z(254 
Z  (  2  5  5 
ZI256 


=6*C*F*R*T*Y 
=B*C*G*R*S*T 
=8*E*G*R*S*T 
=C*£*G*R*S*T 
=8*C*F*R*S#T 
=8*C*G*H*R*Y 
=8*E»G*H*R*Y 
=C*E«G*H*R*Y 
=B*E*F*H*R*Y 
=C*E*F*H*R*Y 
=B»C*6*H#R*S 
=D*E*G*H#R#S 
=8*E*G*H*R*S 
=C*E*G*H#R*S 
*8*EfF*H*R*S 
=C*E*F#H*R*S 
=B*C*G*H*T*Y 
=B*E*G*H*T*Y 
=C*E*G*H*T*Y 
=B*E*F*H*T*Y 
=C*E*F*H*T*Y 
=8*C*G*S*T*Y 
=B*E»G*S*T*Y 
=C*E«G*S*T*Y 
=B*C*F«S*T*Y 
=B*C*G*H*S*T 
=8*E*G«H*S*T 
=C*E*G*H*S*T 
=B*E*F*H*S*T 
=C*E*F*H*S*T 
=8*C*E*F*R*Y 
=B*C*E*F*R*5 
=8*C*E*F*T*Y 
=8*C*E*F*S*T 
=B*C*E*R*S*Y 
=B*E*H*R*S*Y 
=C*E*H*R*S*Y 
=B*C*E*R*T*Y 
=C*B*H*R*T*Y 
=8*E*H*R*T*Y 
=C*E*H*R*T*Y 
=  8  *C*E*R*S*T 
=B*C*H*R*S*T 
=g*F*H*R*S*T 
=C*E»H*R*S*T 
=B*C*E*H*R*Y 
=B*C*E*H*R*S 
=B*C*E*H*T*Y 
=8*C*E*S*T*Y 
=C*B*H*S*T*Y 
=B*E»H*S*T*Y 
=C*E*H*S*>T*Y 
=B*C*E*H»S*T 
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RD=0.0 

DC  3  N=l  .2  56 
C      RD  15  THE  SUM  OF  THE  DENOMINATOR  TERMS 

3  RO=RD+Z(N) 
C      X  IS  THE  TRANSFER  RATIO 
C      AM  IS  THE  TEMPERATURE 

X  «  W  /  R  D 
2  PUNCH  7»AMtX 

END 
C      TYPICAL  INPUT  DATA 


2?3343002 

8  86  2.0 
10000.0 


3293.^ 
10000.0 


845.0 

o.o 


2-500.0 


273.0 
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Table  1 

Voltage  Transfer  Function  for  a  Sample  3-Section  Bridge 

B   =  3,000°  K  Zeros  at   250°  K 

T   =  300°  K  275°  K 

o 

R   =  1.000  ohms  300°  K 


Voltage 

Transfer 

Temperature  °K    Ratio 

225t0C00C 

-.00049 

227.OOOC0 

-.00046 

229.00000 

-.00042 

231 .OO0C0 

-.00039 

233.00000 

-.00035 

235.0OCU0 

-.00031 

237.00000 

-.00u27 

239.COOO0 

-.00023 

241.00000 

-.00018 

243.000CO 

-.00014 

245.00C00 

-.00010 

247.00000 

-.00006 

749.00C00 

-.00002 

251.00000 

.00001 

?53.OdCC0 

.00004 

255.0GCO0 

•000O7 

257. Joi^u 

.00010 

259.0000U 

.00012 

261.00000 

.00013 

263.00C00 

.00013 

265.00000 

.00013 

267.000CO 

.00012 

269.0CC0C 

.00010 

271.0000C 

.00007 

273.00000 

.0  000  4 

275.0CCLL 

J. 00 000 

277.00000 

-.00004 

279.00000 

-.00010 

281 .00000 

-.00015 

2  8  3.^  0 C  0  I 

-.00021 

285.00000 

-.00o26 

2  87.0U000 

-.00030 

289.00000 

-.00033 

291.00C00 

-.00034 

2  9  3.0000  0 

-.00033 

295.0000C 

-.0002  9 

297.C000C 

-.00021 

299.000(  0 

-.00009 

3O1.GUOC0 

•  GO  0'  o  8 

303.00000 

.00032 

Table    1     (continued) 
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3O5.0O0CU 

.00^63 

3U7.U00U0 

.00102 

309.0UC0(j 

.00151 

311  .COOOC! 

.00209 

313.C000C 

.00278 

315.0C00C 

.00360 

317.00CC0 

.00454 

3 19.00CCO 

.00563 

321.00000 

.00686 

'23 .00000 

.00826 

325.00000 

.00982 

69 
APPENDIX  E 
Program  5 

C   PROC-RAM-  TO  COMPUTE  THE  CONTROL  VOL  I  AGE  TO  COMPENSATE  A  CRYSTAL 
DIMMSION  DFR£Q(6o) » CXI 60) .DFRI60) »P( 60 ) . V ( 3 160) iCD(60) 

3  FORMAT  (7F10.1I 

4  FORMAT  (2E10.3) 

5  FOR VAT  (F7.2) 

TMIN  IS  THE  MINIMUM  TEMPERATURE  AND  TMAX  IS  THE  MAXIMUM  TEMPERATURE 
TD  IS  THF  ITERATION  INTERVAL 
FM  AND  FR  ARF  CRYSTAL  CONSTANTS 
FK  IS  A  VARACTOR  CONSTANT 
EP  IS  THE  FREQUENCY  TOLERANCE  IN  PPM 
READ3»TMIN»TMAX,TD.FM,FR»FK>ER 
CO  AND  CF  ARE  CIRCUIT  CAPACITANCES 
REf  i    4,CC,CF 
FUDGE=1Q.**12 
11  T=TM!N 
K  =  l 

DFREO(K)  IS  THE  CRYSTAL  DEVIATION  AT  A  PARTICULAR  TEMPERATURE 
9  READ  5.DFREQ(K) 
VA=-ER 
J  =  I 
DFR(K)  IS  THE  FREQUENCY  PLUS  TCLFRANCE 

6  DFR(K)=-DFREQ(K)+VA 

CX { K )  IS  THE  CRYSTAL  LOAD  CAPACITANCE  FOR  THAT  TEMPERATURE 

CX(K)=-CC-MCC*100tiG0>J.  )/(2.G*FR*<  FM+DFR  (,<))) 

CDCO  IS  THE  VARACTOR  CAPACITANCE  FOR  A.  PARTICULAR  TEMPERATURE 

CDt;<)  =  (CX  (K)-CF  )*2.0 -FUDGE 

V(JiK)  IS  THE  CONTROL  VOLTAGE  FOR  a  PARTICULAR  TEMPERATURE  AND  TOLERANCE 

VtJ»K)=(FK/CD(K) )*»2. 1739130 

J  =  J  +  1 

VA=VA+ER 

IF(FR-VA)  7.6,6 

7  CONTINUE 

PUNCH, T,( V( J,K) , J  =  l  ,3) 
T=T+TD 
<=K.  +  1 

IF(TMAX-T)  1,9,9 
1  CONTINUE 
END 

TYPICAL  INPUT  DATA 
-40.0       70. u        2.0      324. u      208. u      132.0        1.0 
5.00CE-1;  h5.55GE-12 
-A.8r.- 
-?.9r 
-1.10 
.40 
1.6  0 
2.40 
3.41 
4  . 1 C 
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Table  1 

Temperature  °C 

if  If 

-40 

-4.80 

-38 

-2.90 

-36 

-1.10 

-34 

+  .40 

-32 

+  1.60 

-30 

+  2.40 

-78 

+  3.40 

-26 

+  4.10 

-24 

+  4.8  0 

-22 

+  5. 40 

-20 

+  5.80 

-18 

+  6.20 

-16 

+  6.  50 

-14 

+  6.60 

-12 

+  6.70 

-10 

+  6.60 

-  8 

+  6.50 

-  6 

+  6.40 

-  4 

+6.U0 

-  2 

+  5.60 

+  5.20 

+  2 

+  4.80 

+  4 

+  4.20 

+  6 

+  3.60 

+  8 

+  3.O0 

+  10 

+  2.40 

+  12 

+  1.80 

+  14 

+  1  .00 

+  16 

+  .20 

+  18 

-  .60 

+  20 

-1.40 

+  22 

-2.40 

+  24 

-3.30 

+  26 

-4.30 

+  28 

-5.40 

+  30 

-6.40 

+  32 

-7.45 

+  34 

-8.50 

+  36 

-9.40 

+  38 

-10.30 

+  40 

-] 1.10 

+  42 

-12.'.  0 

+  44 

-13.00 

+  46 

-13..  60 

+  48 

-14.40 

+  50 

-15.00 

+  52 

-15.70 

+  54 

-16.30 

+  56 

-16.70 

+  58 

-17.25 

+  60 

-17.80 

+  62 

-18. CO 

+  64 

-18.20 

+  66 

-18.40 

+  68 

-18. 60 

+  70 

-18.70 
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Table  2 

C  C  Temp  °C 

+  1  ppm 

Control 

-1  ppm 

-4. OOOOE+01 

7.5039 

7.5740 

7.644  7 

-3.80UGE+C1 

7.3721 

7.4412 

7.5109 

-3.60U0E+01 

7.24  89 

7.3171 

V. 3859 

-3.400OE+O1 

7.1476 

7.2150 

7.2830 

-3.2C00E+01 

7.1.673 

7.1341 

7.2015 

-3. OOOOE+01 

7.^142 

7.0806 

7.1476 

-2.8000E+01 

6.9482 

7.0142 

7.0806 

-Z.60C CE+O] 

6. 90 24 

6.968" 

7.0341 

-2.4000E+01 

6.8  5  68 

6.9220 

6.9877 

-2.2000E+01 

6.8179 

6.8828 

6.94  8  2 

-2.0000E+0J 

6.792^ 

6.8568 

6.9220 

-1.80GOF+C] 

6.7662 

6.8308 

6.0959 

-1.6000ii+01 

6.74  7'. 

6.8114 

6.6763 

-1.4000E+01 

6.74u6 

6.8049 

6. o698 

-1.200CE+01 

6.7341 

6.7985 

6.0633 

-1.00C0E+01 

6.7406 

6.8049 

6. B698 

-  8  .  0  0  0  C 

6.74  7o 

6.8114 

6.8763 

-6. COCO 

6.7534 

6.8179 

6.8828 

-4.0000 

6.7791 

6.8438 

6.9u89 

-2  .^COfj 

6.8049 

6.8698 

6.9351 

.OCCG 

6.8308 

6.8959 

6.9614 

2.00CO 

6.8568 

6.9220 

6.9877 

4.0000 

6.8959 

6.9614 

7.0274 

6.0000 

6.9351 

7.0C10 

7.0673 

8.0OC0 

6.9746 

7.0407 

7.  1073 

i.ooooe+01 

7.0142 

7.0806 

7. 1476 

1.20U0E+01 

7 . u54o 

7.1207 

7  .  1  6  8  u 

1.4QC0E+01 

7.1073 

7.1745 

7.2421 

1.6000E+01 

7.1610 

7.2286 

7.2966 

1.8000E+01 

7.2150 

7.2830 

7.3514 

2.00GCE+01 

7.2693 

7.3377 

7.4o66 

2.2C0CE+01 

7.3377 

7.4066 

7.4760 

2.400r'E  +  01 

7.3997 

7.4690 

7.5389 

2.6000E+0] 

7.4690 

7.5389 

7.6093 

2.8000E+01 

7.545= 

7.6163 

7.6873 

3. OOOOE+01 

7.6163 

7.6873 

7.7568 

3.2000E+01 

7.6909 

7.7624 

7.o345 

3.4OOGE+01 

7.766u 

7.838  1 

7.9107 

■3.6000E+C1 

7.8309 

7.9034 

7.9  766 

3.8000E+01 

7.8962 

7.9692 

8.0428 

4. OOOOE+01 

7.9546 

fa. 0231 

8.1021 

4.2000E+01 

8.0207 

8.0947 

8.1693 

4.4000E+01 

8.0947 

8.1693 

8.2444 

4.6000E+01 

8.1394 

8.2143 

8.2897 

4.8000E+01 

8.1992 

8.2746 

8.3505 

5.0000E+01 

8.2444 

8.3201 

8.3963 

5.2OCGE+01 

8.29  73 

8.3734 

8  .  4  5  0  0 

5 .4000E+01 

8.3429 

8.4193 

8.4963 

5.6000E+O1 

8.3734 

6.4500 

8.5272 

5 .8000E+01 

8.4154 

8.4924 

8. 5699 

6.0C00E+01 

8.4577 

8.5350 

8.6128 

6.2CO0l:+C1 

8.4731 

8.5505 

8.6284 

6.4C  OCE+01 

8.4885 

8.5660 

8.6441 

6.6000E+01 

8.5040 

8.5816 

8.6598 

6.8000E+0] 

8.5195 

8.5972 

6.6755 

7. OOOOE+01 

8.5272 

8.6050 

8.6833 
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Table  3 

3000°  K    T   ■ 

0 

■  300°  K    R   -  10 
o 

Thermistor- 

.000    R   =  1( 
P 

Temperature 

Resistor 

•k 

Combination 

Thermistor 

1.51C0E+O2 

9.9995E+03 

1 .9294E+08 

1.5300E+02 

9.9993E+03 

1.4881E+08 

1.55COE+02 

9.9991E+03 

1.1555E+08 

1.570CE+02 

9.9989E+03 

9.0298E+07 

1.590OE+Q2 

9.9966E+03 

7.1UO5E+07 

1.6100E+02 

9.9982E+03 

5.6169E+07 

1.6300E+02 

9.9978E+03 

4.4689E+07 

1.65U0E+C2 

9.9972E+C3 

3.S7b3£+07 

1.67G0E+O2 

9.9965E+C3 

2.S756E+C7 

1.6900E+02 

9.9957E+03 

2.3250E+07 

1.7100E+02 

9.9947E+03 

1.B891E+07 

1.730CE+02 

9.9935E+03 

1 .5423E+C7 

1.75O0E+02 

9.992 1E+03 

1.265UE+07 

1.7700E+C2 

9.9904E+03 

1.0423E+07 

1.79C0E+02 

9.9884E+03 

8.6245E+06 

1.81CCE+02 

9.9861E+03 

7.1665E+06 

1.83O0E+P2 

9.9833E+03 

5.9791E+o6 

1.850GE+O2 

'  9.98^1E+03 

5.0^ooE+06 

1.87OCE+02 

9.9763E+03 

4.21o6E+06 

1.8900E+02 

9.9719E+03 

3.5532E+06 

1.9100E+02 

9.9669E+03 

3.0O91E+06 

1.9300E+02 

9.961CE+03 

2.5571E+06 

1.55006+02 

9.9543E+03 

2.18O3E+06 

1.97O0E+02 

9.9467E+03 

1.8650E+C6 

1.9900E+02 

9.9379E+03 

1.6UG3E+06 

2.01C0E+02 

9.9279E+03 

1.3774E+06 

2.03O0E+02 

9.9166E+G3 

1 .1891E+06 

2.O5O0E+02 

9.9038E+03 

1.0294E+C6 

2.0700E+C2 

9.8893E+03 

6.9367E+05 

2.G90GE+C2 

9.8731E+03 

7.7794E+G5 

'2.11O0E+02 

9.8349E+03 

6.789bfc+05 

2.130OE+02 

9.8345E+03 

5.94i3t+05 

2.15O0E+02 

9.8117E+03 

5.2117E+05 

2.1700E+C2 

9.7865E+03 

4.5828E+05 

2.190CE+C2 

9.7584E+03 

4.0392E+05 

2.2100F+02 

9.7274E+03 

3.5683E+05 

2.23C0L+02 

9.6932E+03 

3.1S92E+05 

2.2500E+02 

9.6555E+03 

2.8C/32E+05 

2.27OCE+02 

9.6143E+03 

2.49k5E+C5 

2.29C0E+O2 

9.5691E+C3 

2.22t,8E  +  u5 

2.310GE+02 

9.5198E+03 

1.9826E+05 

2 .33i~0e+G2 

9.4662E+C3 

1 .7  1  jbE  +  Ob 

2.3500E+02 

9.40S1E+03 

1.5894E+05 

2.37OGE+02 

9.3452E+03 

1.4271E+05 

2.391.CE  +  02 

9.2773E+03 

1.2837E+05 

2.41C0E+02 

9.2C42E+03 

1.1567E+05 

10,000 
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2.43G0E+02 
2.4500F+02 
2.470GL+C2 
2.490CE+02 
2.51CGE+02 
2.53GCE+C2 
2.5500E+02 
2. 57CGE+G2 
2.590CE+02 
2.6100E+02 
2.63C0E+C2 
2.650GE+02 
2.6700E+02 
2.6900E+02 
2.710CE+02 
2.730(  E+C2 
2.750UE+02 
2.77G0E+G2 
2.79CCE+02 
2.81C0E+02 
2.83CCE+02 
2.850GE+O2 
2.87C0E+02 
2.8900E+02 
2.910CE+02 
2.9300E+02 
2.950CE+C2 
2.97GCE+02 
2.99G0E+Q2 
3.01GGE+C2 
3.03GCE+C2 
3.G5UCE+C2 
3.07GCE+02 
3.09G0E+02 
3.1 1 OOE+02 
3.1300E+02 
3.15OGE+02 
■3.1700F+02 
3.19OCE+02 
3.21GCE+02 
3.23GCE+02 
3.25GCi.+02 
3.27GCE+G2 
3.29GGE+02 
3.3100E+02 
3.33GGE+C2 
3.35C0E+02 
3.37C0E+G2 
3.39C DE+02 
3.41C0E+02 
3.43C0E+02 
3.4500E+02 
3.470CE+C2 
3.490CE+G2 
3.51LOF+G2 


9.1759E+03 
9.0421E+03 
8.9527E+03 
8.8576E+03 
H.7569E+03 
b.65o3E+03 
8.5380E+03 
6.42t^GE  +  03 
S.2963E+03 
8.1672E+G3 
8.0327E+C3 
7.8931E+03 
7.7486E+03 
7.5995EH  03 
7.4462E+G3 
7.2889E+C3 
7.1281E+03 
6.9642E+03 
6.7976E+G3 
6.6288E+03 
6.4582E+03 
6.2862E+03 
6.1134E+C3 
5.94C2E+C3 
5.7671E+C3 
5.5944E+G3 
5.4227E+C3 
5.2523E+03 
5.0836E+03 
4.917UE+03 
4.7527E+03 
4.5911E+G3 
4.4324E+03 
4.2769E+03 
4.1249E+C3 
3.9763E+03 
3.8315E+03 
3.69G5E+03 
3.5535E+03 
3.42>;4E  +  C3 
3.2914E+C3 
3.166t>E  +  03 
3.C456E  +  0.3 
2.9288E+C3 
2.8160E+03 
2.7072E+03 
2.6023E+03 
2.5013E+03 
2.4C41E+C3 
2.31G6E+03 
2.2207E+03 
2.1343E+C3 
2.0514E+G3 
1.9718E+03 
1.8954E+03 


] .044GE+C5 
9.4395E+G4 
8.5464t+u4 
7.7539E+G4 
7.0441E+04 
6.4u9uE+04 
5.6399E+G4 
5.32o9E+04 
4.8696E+G4 
4.456GE+G4 
4.083GE+G4 
3.7462E+G4 
3.4417E+C4 
3.1656E+04 
2.9157E+G4 
2.6866E+04 
2.4821E+Q4 
2.294lc+G4 
2.12^7E+o4 
1 .9663E+04 
1.8234E+G4 
1.6927E+C4 
1.5730E+04 
1.4632E+04 
1.3624E+04 
1.2699E+G4 
1.1847E+G4 
1.1063E+04 
1.034jL+G4 
9.6732E+03 
9.0373E+03 
6 •486^E+G5 
7.9611E+G3 
7.4732E+G3 
7.02G9E+03 
6.6G12E+C3 
6.2114E+03 
5.8492E+03 
5.5123E+G3 
5.1985E+03 
4.9G63E+03 
4.b337E+C3 
4.37V3E+U3 
4.1416E+G3 
3.9196E+03 
3.7121E+G3 
3.5177E+03 
3.3356E+03 
3.165GE+03 
3.0G-49E+03 
2.8546E+03 
2.7135E+03 
2.58i^8E  +  G3 
2.4361E+G3 
2.3387E+C3 
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3.5300E+02 

1.8221E+03 

2.2201E+03 

3.55OOE+02 

1.7519E+03 

2. 1 24DE+03 

3.57O0E+02 

1.6845E+03 

2.0268E+U3 

3.590CE+02 

1.6ZC0E+03 

1.9331E+03 

3.61OOE+02 

1.5581E+03 

1.8457E+03 

3.630CE+02 

1.4988E+C3 

1.7631E+03 

3.65C0E+02 

1.4420E+03 

1  .685GE+03 

3 .67O0E+02 

1.3876E+03 

1.6112E+03 

3.69GGE+U2 

1.3355E+03 

1.5414E+03 

3.710CE+02 

1.2856E+03 

1.4753E+C3 

3.7300E+02 

1.2378E+03 

1.4127E+03 

3.7&'JGE  +  U2 

1.1920E+03 

1.3534E+03 

3.7700E+02 

1.1482E+03 

1.2971E+03 

3.79C0E+O2 

1.1C62E+03 

1.2436E  +  CJ3 

3.8100E+02 

1.0660E+Q3 

1.1932E+03 

3.830CE+02 

1.0274E+03 

1.1451E+03 

3.8500E+02 

9.9053E+02 

1.0994E+03 

3.870CE+C2 

9.5517E+02 

1.056GE+03 

3.89OOE+02 

9.2129E+02 

1.014&E+03 

3.91GCE+02 

8.8882E+G2 

9.7552E+D2 

3.930CE+02 

8.5770E+02 

9.3817E+02 

3.95GGE+G2 

8.2787E+02 

9.026GE+02 

3.97UOE+02 

7.9928E+02 

8.6872L+G^ 

3.9900E+02 

7.7186E+02 

8.3643E+02 

4.01G0E+02 

7.4557E+02 

8.0564E+C2 

4.0300E+02 

7.2036E+02 

7.7628E+02 

4.0500E+02 

6.9617E+02 

7.4826E+02 

4.0700E+02 

6.7296E+02 

7.2151E+02 

4.O9O0E+C2 

6.5068E+02 

6.9597E+02 

4.11GCF+C2 

6.2930E+02 

6.71t6E+u2 

4.13Q0E+02 

6.0878E+02 

6.4824E+G2 

4.15U0E+02 

5.89u7E+02 

6.2594E+02 

4. 17G0E+G2 

5.7G14E+C2 

6.C461E+G2 

4.19UGE+C2 

5.51V6E+02 

b.642uE+G2 

4.21C0E+02 

5.3446E+02 

3.6466E+G2 

4.230CE+C2 

5.1769E+02 

5.4596E+02 

4.250C  i>02 

5.0155E+02 

5.28u4£+02 

4.2700E+02 

4.8603E+02 

5.1086E+02 

4.?9rCE+r2 

4.711 1E+02 

4.9440E+02 

4.3100E+02 

4.5675E+02 

4.7361E+02 

4.33C OE+02 

4.4294F+C2 

4.6347E+02 

4.350GE+02 

4.2965E+Q2 

4.4894E+02 

4.3700E+02 

4.1636E+C2 

4.3499E+02 

4.39J0E+02 

4.0454E+02 

4.216GE+C2 

4.4100E+02 

3.9266E+02 

4.0b73E+02 

4.43C0E+02 

3.8126E+02 

3.9637E+02 

4.45i-GE  +  u2 

3.7t,25E  +  C2 

3.8449E+02 

4.47GCE+02 

3.5965E+02 

3.73G6E+C2 

4.4900E+02 

3.4342E+G2 

3.62U&E+U2 

B  =  2,500"  K 


Temperature 

•K 

1.51lC_  +  02 
1.53COE+02 

1.55O0F+02 
1.570CE+02 
1.5900E+02 
1.61O0E+02 
1.6300E+C2 
1.65OCE+02 
1.670CE+C2 
1.6900E+02 
1.71OGE+02 
1.73C0E+02 
1.7500E+02 
1.7700E+02 
] .790CE+02 
1.8100E+02 
1.8300E+02 
1.8500E+02 
1.87CCE+02 
1.89CCE+02 
1.910GE+C2 
1.930GE+02 
1.950CE+02 
1.97CCE+02 
1.99GOE+02 
2.01C0E+02 
2.030OE+02 
2.O5OOE+02 
2.070OE+02 
2.C900E+02 
2.110CE+02 
2.13C0E+02 
2.15UCE+02 
2.17UCE+C2 
2.19O0E+02 
2.21O0E+02 
2.23CCE+02 
2.2500E+02 
2.27COE+02 
2.290CL+02 
2.31C0E+O2 
2.33GOE+02 
2.35UOE+C2 
2.370CE+02 
2.39LCE+C2 
2.410CE+Q2 
2.43O0E+O2 
2.45U0E+C2 
2.47O0E+O2 
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Table 

4 

3°  K    R   = 
o 

10,000    R   -  10,000 
P 

Thermistor- 

Resistor 

Combinat  ion 

Thermistor 

9.9939E+03 

9.9924E+03 

1.634GE+07 
1.3160E+07 

9.9906E+03 

1.0658E+C7 

9.9885E+03 

8.6788E+G6 

9.9859E+03 

7.1G35E+06 

9.9829E+03 

5.&432E+G6 

9.9793E+C3 

4.829bE+06 

9.97blE+03 

4.01u2E+u6 

9.97L.2E  +  03 

3.3447E+J6 

9.9644E+03 

2.8ui7E+u6 

9.9577E+03 

2 .3566E+ j6 

9.95O0E+03 

1.99G1E+C6 

9.9411E+03 

1.6871E+06 

9.9308E+C3 

1 .4356E+06 

9.9I91E+03 

1.2261E+06 

9.9C57E+03 

1 .05G7E+C6 

9.89u5E+03 

9.03bGE+u5 

9  .  8  7  3  3  E  +  0  3 

7.7945E+05 

9.8539E+03 

6.7436E+0S 

9.8321E+03 

5.85S7E+05 

9.8L76E+03 

5.09B3E+C15 

9.78u3E+03 

4.4516E+05 

9.7499E+03 

3.8978E+05 

9.7161E+03 

3.4221E+G5 

9.6787E+03 

3.0123E+05 

9.6375E+03 

2.65B4E+05 

9.5921E+03 

2.3518E+05 

9.5424E+03 

2.0855E+05 

9.4882E+03 

1.8537E+05 

9.4290E+03 

1.6514E+05 

9.3648E+03 

1.4744E+G5 

9.2954E+03 

1 .3192E+05 

9.2204E+03 

1.1827E+05 

9.1398E+03 

1 .06<^5E  +  u5 

9.C534E+03 

9.5641E+04 

8.9611E+03 

8.62b4E+04 

-' 

8.8628E+03 

7.7932E+G4 

8.7584E+03 

7.0541E+O4 

8.648r:E  +  C3 

6.3962E+04 

8.5315E+G3 

5.8u97E+04 

8.4C91E+03 

5.2857E+04 

3.28O8E+03 

4.3166E+C4 

8.1469E+03 
8.0G7bE+03 
7.8630E+03 
7.713bE+03 
7.bS94E+G3 
7.4C12E+03 
7.2391E+03 


4.3964E+04 
4.C1S9E+04 
3.679<4E  +  u4 
3.3734E+C4 
3.C974E+G4 
2.O479E+04 
2.622uE+04 
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Table  4  (continued) 

2.4900E+02 

7.0737E+C3 

2.4173F+04 

2.51OOE+02 

6.9054E+03 

2.2314F+04 

2.5300E+02 

6.7347E+03 

2.Q625E+04 

2.55OOE+02 

6.5620E+03 

1.9087E+04 

2. 5700E+02 

6.3S79E+03 

1 .7665E+G4 

2.5900E+02 

6.2128E+03 

1 .64G6E+04 

2.61L.CL+C2 

6.0373E+03 

1.5235E+04 

2.63OUE+02 

5.8618E+03 

1.41bbt+G4 

2.6500E+02 

5.6868E+03 

1 .31b4E+G4 

2.6700E+C2 

5.5127E+03 

1.2285E+04 

2.6900E+02 

5.3399E+03 

1.14i>9E  +  G4 

2.71 tOF +02 

5.1689E+C3 

1.0699E+G4 

2.7300E+02 

5.00OCE+O3 

1  .0O00E+04 

2.7500E+02 

4.8336F+03 

9.3557E+03 

2.7700E+O2 

4.6699F+C3 

8.7613E+03 

2.7900E+02 

4.5092E+03 

8.2124E+C3 

2.81O0E+02 

4.3519E+03 

7.7O50E+03 

2.8301.F  +  C2 

4.198GE+03 

7.23&5E+G3 

2.8500E+02 

4.0476E+C3 

6.8ou6t+03 

2.87O0E+02 

3.9G14E+03 

6.3973E+03 

2.890OE+02 

3.7590E+03 

6.0231E+u3 

2.9100E+02 

3.62 G6E+C3 

5.67S4E+Q3 

2.9300E+02 

3.4863E+03 

5.3522E+03 

2.9500E+02 

3.3561E+03 

5.0513E+U3 

2.9700E+02 

3.2300E+03 

4.7711E+03 

2.990CE+02 

3. 1082E+03 

4.S'.99E  +  C3 

3.0100E+02 

2.99G4E+03 

4.2662E+03 

3.O3O0E+O2 

2.876SE+C3 

4.0386E+03 

3.C500E+O2 

2.7672E+03 

3.8259E+03 

3.070CE+C2 

2.6616E+C3 

3.627GE+03 

3.U900E+02 

2.5b99E+03 

3.44O6E+03 

3. 1100E+02 

2.4621E+03 

3.2663E+03 

3.13OCE+02 

2.3680E+03 

3.1^2bE+G3 

3.1500E+02 

2.2776E+C3 

2.9493E+03 

3.1700E+02 

2.1907E+03 

2.8053E+03 

3.1900E+T2 

2.1073E+03 

2.6700E+03 

3.2100E+02 

2.0273E+03 

2.5427E+03 

3.23O0E+02 

1.9504E+03 

2.423uE+G3 

3.250OF+O2 

1.S767E+C3 

2.3103E+03 

. 

3.27J0E+02 

1.8061E+03 

2.2G41E+03 

3.29C0E+C2 

1.7383E+03 

2.1^40E+03 

3.310CE+C2 

1.6733E+03 

2.0O96E+03 

3.3300E+02 

1.6111E+03 

1.9205E+03 

3.3500E+02 

1.5514E+C3 

1.8363E+Q3 

3.370GF+C2 

1.4943E+03 

1 .736oE+03 

3.39C0i.  +  C2 

1.4395E+03 

1.6815E+03 

3.410C  E+02 

1.387GE+03 

1.6104E+03 

3.430CE+02 

1.3367E+03 

1.543CE+03 

3.4500E+02 

1.2885E+03 

1.4791E+03 
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3.4700E+02 

1.2424E+03 

1.4186E+03 

3.4900E+02 

1.1981E+03 

] .3612E+03 

3.51G0E+G2 

1.1558E+03 

1.3u68E+03 

3.5300E+02 

1.1152E+03 

1.2551E+03 

3.5500E+C2 

1.0762E+03 

1.2060E+03 

3.57O0E+02  ■ 

1.0388E+03 

1.1564E+03 

3.5900E+02 

1.0032E+G3 

1.1150E+03 

3.61U0I.+G2 

9.6887E+02 

1.0728E+03 

3.6300E+02 

9.3599E+02 

1.0326E+03 

3.65OOE+02 

9.0447E+02 

9.9441E+02 

3.6700E+02 

3.7423E+G2 

9.5797E+G2 

3.6900E+02 

8.4521E+02 

9.2325E+G2 

3.73  00E+02 

8.1738E+02 

8.9014E+02 

3.7300E+02 

7.9067E+02 

8.5855E+02 

3.75GGB+02 

7.65G3E-»C2 

8.284GE+02 

3.77G0E+G2 

7.4G41E+02 

7.9962E+02 

3.79GGE+G2 

7.1677E+02 

7.7212E+C2 

3.81O0E+02 

6.9407E+02. 

7.4584E+02 

3.83LGE+G2 

6.7227E+G2 

7.2U72E+G2 

3.8500E+02 

6.5131E+02 

6.9669E+C2 

3.87OOE+02 

6.3117E+02 

6.737GE+02 

3.89GOF+02 

6.1182E+02 

6.5169E+C2 

3.9100E+02 

5.9321E+02 

6.3O61E+02 

3 .9300F+02 

5.7531E+02 

6.1042E+C2 

3.9500E+02 

-  5.5809E+02 

5.9108E+C2 

3.9700E+02 

5.4152E+02 

5.7253E+02 

3.9900E+02 

5.2558E+02 

5.5474E+02 

4.0100E+02 

5.1G24E+02 

5.3767E+02 

4.C3G0E+C2 

4.9546E+02 

5.2129E+02 

4.05U0E+02 

4.8123E+02 

5.0bb6t+G2 

4.G7GGE+02 

4.6753E+02 

4.9U46E+02 

4.G9GGE+G2 

4.5432E+02 

4.7594E+02 

4.1100E+G2 

4.4159E+02 

4.6200E+02 

4.13GGE+02 

4.2933E+02 

4.4859E+02 

4.1500E+02 

4.] 75GE+02 

4.3569E+02 

4.  17006+02 

4.06G9E+02 

4.2328E+02 

4.1900E+02 

3.95u9E+G2 

4.1134E+02 

4.2100E+02 

3.8447E+02 

3.99o4E+G2 

4.23GGE+02 

3.7422E+02 

3.8877E+G2 

4.2500E+02 

3.6433E+02 

3.7811E+G2 

4.2700E+02 

3.5478E+02 

3.67&3E+G2 

4.2900E+02 

3.4556E+02 

3.5793E+02 

4.3100E+G2 

3.3665E+Q2 

3.48;>6E  +  G2 

4.33C0E+02 

3.2SG4E+02 

3.3917E+G2 

4.350CE+02 

3.1973E+02 

3.3U29E+02 

4.3700F+C2 

3.1168E+C2 

3.2171E+02 

4.3900L+02 

3.0391E+02 

3.1344E+G2 

4.4100E+02 

2.9639E+C2 

3.0544E+02 

4.43C0E+02 

2.8912E+02 

2.9773E+02 

4.4500E+02 

2.82C8E+02 

2.9'-27£  +  G2 

4.47GCE+G2 

2.7527E+02 

2.83G6E+02 

4.49CCF+G2 

2.6868E+02 

2.761GE+02 
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Table  5 

3-Section 

Bridge  Transfer  Fu 

Temperature  °K 

Trans.  Function 

233.00C0 0 

.00220 

235.CO00G 

.00267 

237.00000 

.00314 

?3S.0uu^u 

.00360 

241 .uutLo 

.00405 

243.0000U 

.00449 

245.00000 

.00492 

247.0OCG0 

.005 32 

?49.ooooc 

.00570 

;;51  .00000 

.00605 

253.00000 

.00636 

255. 000  0  U 

.00663 

257.00000 

.00686 

259.OG0G0 

. G07G4 

261.000U0 

.00717 

263.GGGOO 

.00724 

265.OUO00 

.00725 

267.00000 

.0072  1 

269.00000 

.00710 

271.00COC 

.00693 

273.0001  0 

.00669 

275.0000C 

.00639 

277.GOCL0 

.00603 

279.00000 

.00561 

281.00000- 

.00513 

283.00000 

.00459 

285.00000 

.00400 

287.OU0OG 

.00336 

2  89 • uuG^u 

.00267 

291.00000 

.00195 

293.00000 

.00119 

295.00000 

.00040 

297.00000 

-.00041 

299.00000 

-.00124 

301.00000 

-.00208 

303.00000 

-.00293 

1 0  5  .  0  0  C  i/  J 

-.00378 

307.0UOOO 

-.00462 

3u9.^GoG0 

-.00545 

311.00000 

-.00625 

313.0o000 

-.00704 

3  15.00000 

-.00776 

917.000CC 

-.00649 

319.00000 

-.00916 

?21.0000 0 

-.00977 

323.00000 

-.01033 

9  25.0 0 C 0 0 

-.0  1084 

^27.00000 

-.01127 

329.00000 

-.01164 

331.00CC0 

-.01193 

333.ouGi.G 

-.01215 

335.00000 

-.01228 

337.000^0 

-.01233 

339.CC000 

-.01230 

341.00000 

-.0121 7 

343.00 C 0 0 

-."1195 
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Table  6 
Bridge  Output  Amplified  and  Shifted 


Temperature  °K 
2.33U0E+02 
2.350CE+02 
2.37GCE+02 
2.39vjGE  +  02 
2.410GE+02 
2.43GOE+02 
2.45OCE+02 
2.470GE+02 
2.4900E+02 
2.510CE+02 
2.5300E+C2 
2.5500E+G2 
2.570GE+02 
2.59GGE+G2 
2.61UGE+02 
2.63i,0E  +  G2 
2.65O0E+O2 
2.67GCE+02 
2.69O0E+C2 
2.710CE+02 
2.7300E+02 
2.7500E+02 
2.7700E+02 
2.790UE+G2 
2.8100E+62 
2.83U0E+02 
2.85uOE+02 
2.8700E+02 
2.89U0E+02 
2.9100E+02 
2.930CE+C2 
2.9500E+02 
2.970CE+02 
2.990CE+C2 
3.Q10CE+02 
3.03L0E+C2 
3.05U0E+02 
3.07GCE+02 
3  .09OGE+C2 
3.  lluGE+02 
3.1300E+02 
3.15O0E+O2 
3.17O0E+O2 
3.19C0E+02 
3.21C0E+02 
3.23'jCE  +  C2 
3.25OOE+02 
3.27C0E+02 
3.29U0E+C2 
3.31GGE+02 
3.  33^0E  +  l.2 
3.  35GGE+G2 
3.37GCE+02 
3.39C0E+02 
3.410OE+O2 
3.4300E+02 


Control  Voltage 
7.276* 
7.2347 
7.1930 
7.152  2 
7.1123 
7.0733 
7.0352 
6.9997 
6.966  0 
6.9349 
6.9075 
6.8835 
6.8631 
6.8471 
6.8356 
6.8294 
6.8285 
6.832] 
6.8418 
6.8569 
6.8782 
6.9048 
6.9367 
6.9740 
7.C165 
7.0644 
7.1168 
7.1735 
7.2  34  7 
7.2986 
7.3660 
7.4  360 
7.5079 
7.5815 
7.656  0 
7.7314 
7.8067 
7.8812 
7.954b 
e.0256 
b.0959 
8.1615 
8,22*5 
6.2839 
8.338'i 
6.3876 
8.4329 
8.4-710 
b.5038 
6.5295 
8.5490 
8.5606 
8.5650 
8.5624 
8.5506 
8.5313 
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Certain  temperature  compensation  problems  require  the 
generation  of  a  temperature  compensated  voltage.   Compensating 
a  quartz  crystal  resonator  for  frequency  variations  is  one  such 
problem.   The  frequency  can  be  controlled  by  loading  the  crystal 
with  a  variable  capacitance.   This  capacitance  can  be  supplied 
by  a  varactor  biased  by  a  temperature  variable  control  voltage. 

The  low  thermistor  cost  and  the  large  variety  of  thermistor 
types  available  make  it  an  excellent  compensation  element.   The 
power  dissipation  constant  of  a  thermistor  must  be  chosen  to 
match  similar  thermal  constants  of  the  compensated  device. 
Thermistor  resistance  is  given  by 


i(l/T  -  1/T  ) 


R  =  R   e 


A  power  series  representation  of  this  equation  was  found  to  be 
impractical.   A  fixed  resistor  was  placed  in  parallel  with  the 
thermistor.   This  combination  could  then  be  approximated  by  a 
linear  expression  over  a  particular  temperature  range.   This 
linear  expression  was  used  in  the  analysis  and  design  of  possi- 
ble compensation  circuits.   Two  networks,  the  cascaded  bridge 
and  the  resistance  ladder  were  examined  as  possible  compensation 
networks.   The  ladder  network  did  not  have  sufficient  flexi- 
bility.  However,  the  cascaded  bridge  network  was  developed  and 
demonstrated  to  be  a  versatile  method  of  generating  temperature 
compensated  voltages.   The  transfer  function  of  the  bridge  net- 
work was  found  by  the  use  of  topological  formulas.   Digital  com- 
puter programs  were  used  as  an  aid  to  the  solution  of  these 


formulas  and  other  calculations.   Voltage  shifting  and  amplifi- 
cation techniques  were  used  as  an  aid  to  the  synthesis  of 
temperature  compensated  voltages.   These  techniques  along  with 
the  cascaded  bridge  network  were  applied  to  the  problem  of  com- 
pensating a  crystal  resonator  for  temperature  induced  frequency 
variations.   The  results  of  theoretical  computations  using  this 
network  show  that  frequency  deviation  can  be  reduced  to  +1  ppm 
or  less  over  a  temperature  range  from  -35°  to  +70°  centigrade. 


